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Abstract 

Polarimetric-hyperspectral  remote  sensing  is  a  promising  field  that  brings  two 
traditionally  independent  modalities  together  to  enhance  scene  characterization 
capabilities.  The  Telops  polarimetric-hyperspectral  imager  (P-HSI)  combines  these 
technologies  and  provides  a  combined  imaging  and  spectral  capability  that  is  considered 
state-of-the  art.  The  Defense  Threat  Reduction  Agency  (DTRA)  funded  research  at  AFIT 
to  leverage  this  capability  to  better  inform  radiation  transport  models  that  will  allow  for 
more  effective  location  of  ionizing  radiation  sources  within  the  scene. 

The  desire  for  highly  accurate  data  requires  careful  calibration  of  the  instrument. 
This  thesis  develops  a  mathematical  calibration  framework  that  links  standard  Fourier 
transform  spectrometer  (FTS)  calibration  with  standard  polarimetric  calibration  in  a 
simple,  straightforward  manner.  Building  this  framework  was  key  for  two  reasons.  First, 
it  was  essential  for  ensuring  the  calibration  methodology  was  sound  and  it  provided  a 
framework  for  understanding  the  influence  of  various  instrument  parameters  (both  ideal 
and  non-ideal)  on  ultimate  calibration  performance.  Second,  this  development  makes  it 
easy  for  both  FTS  and  polarimetric  experts  to  understand,  and  paves  the  way  for  use  of 
this  method  for  similar  instruments  in  the  future. 

This  framework  is  then  utilized  to  quantify  the  non-idealities  of  the  system  and  to 
characterize  the  performance  of  the  spectro-polarimetric  calibration.  The  calibration  was 
first  performed  under  the  assumption  of  an  ideal  polarizer,  and  then  an  additional 
correction  was  made  for  the  non-ideal  nature  of  the  polarizer.  The  results  showed  that  the 
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calibration  was  accurate  to  within  1.18%  and  accounting  for  the  non-ideal  nature  of  the 
polarizer  provided  an  additional  0.56%  correction  in  the  measured  degree  of  linear 
polarization  (DoLP).  It  was  determined  that  the  residual  error  in  the  measured  DoLP  was 
primarily  due  to  a  bias  in  the  data  caused  by  the  thermal  instability  of  one  of  the 
blackbody  sources  and  not  because  of  the  calibration  itself.  Other  sources  of  residual 
error  included  noise,  detector  non-linearity,  and  other  experimental  setup  limitations. 

The  band-averaged  noise  equivalent  spectral  DoLP  (NESDoLP)  was  determined  to  be 
0.39%  at  16  cm"1  resolution  with  a  200  //s  integration  time  and  flood  illumination  from  a 
55  °C  blackbody.  This  represents  the  limit  in  the  accuracy  of  the  calibration  that  should 
be  achieved  after  the  bias  is  removed.  The  band-averaged  noise  equivalent  spectral 
radiance  (NESR)  values  were  30.6  nW/(cm2  sr  cm"1),  33.6  nW/(cm2  sr  cm'1), 

33.6  nW/(cm2  sr  cm"1),  and  36.9  nW/(cm2  sr  cm'1)  for  the  0°,  45°,  135°,  and  90°  polarizer 
measurements  respectively.  This  represents  the  limit  in  the  accuracy  of  the  radiometric 
calibration.  The  diattenuation  of  the  hyperspectral  imaging  system,  excluding  the 
polarizer,  was  determined  to  be  0.43  which  corresponds  to  an  extinction  ratio  of  2.51. 
This  is  small  compared  with  the  polarizer's  extinction  ratio  of  400;  however,  it  must  be 
calibrated  out  to  achieve  the  most  accurate  data  from  the  instrument.  This  thesis 
demonstrated  the  capability  to  perform  radiometric,  spectral,  and  polarimetric  calibration 
in  order  to  achieve  highly  accurate  data  from  the  Telops  instrument.  This  will  enable 
more  accurate  scene  characterization,  which  will  better  inform  the  radiation  models 
developed  under  DTRA  sponsorship. 
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Polarimetric  Calibration  and  Characterization  of  the  Telops  Field  Portable 
Polarimetric-Hyperspectral  Imager  in  the  Long  Wave  Infrared 


I.  Introduction 


Overview 

Polarimetric-hyperspectral  remote  sensing  is  a  promising  field  that  brings  two 
traditionally  independent  modalities  together  to  enhance  scene  characterization 
capabilities.  The  Telops  polarimetric-hyperspectral  imager  (P-HSI)  combines  these 
technologies  and  provides  a  combined  imaging  and  spectral  capability  that  is  considered 
state-of-the  art.  The  Defense  Threat  Reduction  Agency  (DTRA)  funded  research  at  AFIT 
to  leverage  this  scene  characterization  capability  to  provide  better  information  to 
radiation  transport  models  that  will  allow  for  more  effective  location  of  radiation  sources 
within  a  region  of  interest.  As  a  baseline,  the  models  require  as  much  information  about 
the  region  of  interest  as  possible,  including  the  materials  in  the  scene  and  the  orientation 
of  objects  within  the  scene.  Combining  the  hyperspectral  and  spectro-polarimetric  data 
sets  in  a  system  with  high  spatial  and  spectral  resolution  provides  a  more  robust  scene 
characterization  capability  that  will  be  utilized  in  the  radiation  models. 

To  support  the  objectives  of  the  DTRA  effort,  there  is  a  requirement  for  highly 
accurate  radiometric,  polarimetric,  and  spectral  data  on  a  pixel-by-pixel  basis.  The 
complex  nature  of  the  Telops  instrument  combined  with  working  in  the  thermal  IR 
waveband  makes  achieving  this  accuracy  a  challenge.  Because  of  this,  great  care  must  be 
taken  to  develop  a  sound  calibration  methodology.  This  thesis  develops  and 
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demonstrates  the  ability  to  perform  calibration  in  all  three  domains.  That  is,  radiometric 
calibration,  off-axis  spectral  calibration,  and  polarimetric  calibration  are  now  successfully 
implemented.  This  means  that  highly  accurate  radiometric,  spectral,  and  polarimetric 
data  can  be  generated  with  the  Telops  instrument.  This  will  enable  more  accurate  scene 
characterization,  which  will  better  inform  the  radiation  models  developed  under  DTRA 
sponsorship. 

Additionally,  this  thesis  develops  a  mathematical  calibration  framework  that  links 
standard  Fourier  transform  spectrometer  (FTS)  calibration  with  standard  polarimetric 
calibration  in  a  simple,  straightforward  manner.  Building  this  framework  was  key  for 
two  reasons.  First,  it  was  essential  for  ensuring  the  calibration  methodology  was  sound 
and  it  provided  a  framework  for  understanding  the  influence  of  various  instrument 
parameters  (both  ideal  and  non-ideal)  on  ultimate  calibration  performance.  Second,  this 
development  makes  it  easy  for  both  FTS  and  polarimetric  experts  to  understand,  and 
paves  the  way  for  use  of  this  method  for  similar  instruments  in  the  future. 

Organization  of  Thesis 

Chapter  II  of  this  thesis  starts  off  by  providing  an  overview  of  the  basic  theory 
and  concepts  required  to  understand  the  technology  that  comprises  the  Telops  instrument 
as  well  as  the  methodology  described  in  later  chapters.  Next,  Chapter  III  covers  in  some 
detail  the  radiometric  calibration  and  off-axis  spectral  calibration  that  are  required 
whether  the  polarizer  is  there  or  not,  and  provides  results  which  demonstrate  the  effect 
that  these  calibrations  have  on  the  quality  of  the  data.  This  chapter  is  a  useful  primer  for 
the  full  spectro-polarimetric  calibration  that  is  developed  in  Chapter  IV.  Chapter  IV  is 


2 


considered  the  primary  chapter  of  this  thesis.  It  develops  the  full  spectro-polarimetric 
calibration  of  the  instrument,  providing  a  mathematical  framework  that  outlines  the 
calibration  methodology.  This  methodology  is  then  used  to  quantify  the  non-idealities  of 
the  system  and  to  characterize  the  performance  of  the  spectro-polarimetric  calibration  on 
experimental  data.  Additionally,  the  noise  equivalent  spectral  radiance  (NESR)  and  noise 
equivalent  spectral  degree  of  linear  polarization  (NESDoLP)  are  determined  for  the 
dataset,  which  indicates  the  fundamental  limit  of  performance  achieved  with  the  system. 
In  Chapter  V,  the  polarimetric  system  model  is  then  used  to  simulate  the  effect  that 
changing  various  system  parameters  has  on  the  measured  data  and  is  used  to  estimate  the 
diattenuation  and  extinction  ratio  of  the  hyperspectral  imager  itself  (not  including  the 
polarizer).  Finally,  Chapter  V  provides  the  conclusions  and  future  recommendations. 
These  recommendations  provide  useful  ways  to  improve  the  accuracy  of  the  calibration. 

Benefits  of  Study 

Beyond  the  benefits  mentioned  above,  this  research  provides  broad  benefits  to  the 
Department  of  Defense  (DoD)  and  intelligence  community  (IC)  in  general.  Figure  1 
outlines  some  of  the  benefits  of  each  technology.  Target  detection/material 
identification,  anomaly  detection,  change  detection,  and  buried  improvised  explosive 
device  (IED)  detection  are  all  capabilities  of  high  interest  to  the  DoD  and  IC. 
Furthermore,  the  ability  to  do  so  in  day  or  night  when  the  target  is  camouflaged  or  hidden 
within  a  natural  background  and  potentially  occupying  very  few  pixels  within  the  scene  is 
a  high  priority.  Hyperspectral,  polarimetric,  and  long-wave  infrared  (LWIR)  imaging  all 
provide  a  unique  capability  to  discriminate  a  target  from  background  scenery. 
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Hyperspectral  imaging  is  useful  for  characterization  of  materials  within  the  scene.  This 
allows  the  possibility  of  discriminating  a  target  based  on  its  spectral  characteristics. 
Polarimetric  imaging  is  useful  for  detecting  manmade  objects  in  a  natural  background 
since  manmade  objects  tend  to  have  different  polarimetric  characteristics  than  natural 
objects.  It  has  also  been  shown  to  provide  a  capability  for  disturbed  earth  detection  and 
therefore  buried  mine/IED  detection.  LWIR  imaging  is  useful  for  detecting  objects  in 
both  day  and  night  that  are  uncharacteristically  warm  or  cool  with  respect  to  the 
surrounding  environment  whether  the  object  is  concealed/camouflaged  or  not.  The 
combination  of  these  technologies  provides  an  enhanced  dataset  that  contains  the 
traditional  hyperspectral  signatures  as  well  as  spectrally  resolved  polarimetric  signatures. 
This  combination  enables  a  more  robust  knowledge  of  the  complex  index  of  refraction, 
which  is  a  fundamental  material  property  that  governs  light  propagation,  reflection,  and 
absorption  within  the  material.  This  results  in  an  improved  ability  to  characterize  the 
scene  and  provide  the  capabilities  listed  in  Figure  1 . 
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Figure  1.  Flow  chart  outlining  some  of  the  benefits  of  each  technology  and  the  enhanced  capability 

enabled  by  the  combination  of  the  technologies. 
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II.  Background  and  Theory 


In  order  to  understand  how  the  P-HSI  instrument  works,  it  is  useful  to  understand 
how  each  component  technology  works.  This  section  is  a  review  of  the  basic  concepts 
associated  with  each  major  technology  used  in  the  Telops  imager. 

Introduction  to  Hyperspectral  Imaging 

Hyperspectral  imaging  has  a  wide  range  of  commercial  applications  including 
agriculture,  environmental  monitoring,  and  mineral  exploration.  The  DoD/IC  also  use 
hyperspectral  imaging  in  a  variety  of  ways  including  characterization  of  effluents  coming 
from  facilities,  background  classification,  target  detection,  anomaly  detection,  and  change 
detection  using  overhead  airborne  or  space-based  platforms.  Hyperspectral  imaging 
involves  a  combination  of  imaging  and  spectroscopy  technologies.  Imagery  technology 
allows  the  measurement  of  light  intensity  over  a  given  spatial  domain.  Spectroscopic 
technology  allows  the  measurement  of  light  intensity  as  a  function  of  wavelength. 
Combining  the  two  technologies  enables  a  measurement  of  light  intensity  with  respect  to 
both  the  spatial  dimension  and  the  spectral  dimension.  Therefore  each  pixel  of  the 
camera  contains  not  just  intensity,  but  intensity  as  a  function  of  wavelength.  This 
concept  is  summarized  in  Figure  2. 
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Figure  2.  Illustration  of  the  hyperspectral  imaging  data  cube.  Each  layer  represents  a  picture  of  the 
scene  for  a  given  spectral  band.  For  a  typical  color  image  this  would  contain  three  layers,  one  for 
red,  green,  and  blue.  For  a  hyperspectral  image,  this  contains  hundreds  of  layers.  Each  pixel 
contains  a  spectrum  that  can  be  used  to  identify  the  material  present  in  the  pixel.  (Image  adapted 

with  permission  from  Ref  [1]). 


As  Figure  2  demonstrates,  the  result  of  a  collection  of  hyperspectral  images  is 
typically  represented  as  a  data  cube.  This  data  cube  contains  the  spatial  information  in 
the  x-y  plane  and  the  spectral  information  in  the  z-direction.  Since  each  pixel  contains 
hyperspectral  information,  it  is  possible  to  identify  a  target  of  interest  even  when  that 
target  only  occupies  one  pixel.  With  the  use  of  proper  post  processing,  this  provides  a 
great  target  detection  capability  that  is  not  achieved  with  standard  color  imagery.  One  of 
the  challenges  with  this  is  the  large  amount  of  data  produced  by  high  resolution 
hyperspectral  imagery.  Modern  technology  has  enabled  the  ability  to  process  such  large 
data  sets  efficiently  and  thus  enabled  its  use  in  military  remote  sensing  applications 
where  in  the  past  it  was  not  feasible. 
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The  capability  to  do  material  identification  with  hyperspectral  remote  sensing  is 
made  possible  by  the  quantum  mechanical  nature  of  the  interaction  of  light  with  matter. 
All  matter  consists  of  atoms  and  molecules  and  it  is  the  atomic  and  molecular 
spectroscopy  derived  from  quantum  mechanics  that  give  the  matter  its  optical  and 
spectral  properties.  These  properties  result  in  spectral  features  that  are  observed  during 
hyperspectral  remote  sensing  that  can  potentially  act  as  a  finger  print  for  the  material 
being  observed.  A  physical  parameter  that  defines  the  interaction  of  radiation  with  matter 
is  the  complex  index  of  refraction,  given  by 

n  =  n  +  i/c  (1) 

where  the  real  part,  n ,  determines  the  phase  propagation  in  the  material  and  the 
imaginary  part,  k  ,  determines  the  absorptive  properties  of  the  material.  It  is  the  spectral 
dependence  of  this  intrinsic  material  property  that  determines  the  spectral  signatures 
observed  in  hyperspectral  remote  sensing.  More  specifically,  it  is  k  that  governs  the 
spectral  absorption  and  emission  features  of  the  material,  and  therefore,  its  spectral 
signature  observed  in  traditional  hyperspectral  remote  sensing.  With  the  addition  of 
spectrally  resolved  polarimetric  data,  more  information  about  the  complex  index  of 
refraction  can  be  determined  as  well  through  the  Fresnel  equations  and  the  proper 
polarimetric  bi-directional  reflectivity  distribution  (pBRDF)  modeling. 

Fourier  Transform  Spectrometer 

The  instrument  that  captures  the  spectral  information  for  the  image  is  the  Fourier 
Transform  Spectrometer  (FTS).  The  FTS  in  the  Telops  is  a  Michelson  interferometer  that 
uses  corner  cubes  rather  than  parallel  mirrors.  In  the  Michelson  interferometer,  a  beam 
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splitter  is  used  to  divide  the  light  into  two  beams.  After  the  two  beams  have  each 
reflected  off  of  a  different  mirror,  the  two  beams  are  recombined  by  the  same  beam 
splitter  and  sent  to  the  detector.  One  of  the  mirrors  is  fixed  and  one  is  movable.  This 
allows  the  instrument  to  change  the  optical  path  difference  (OPD)  between  the  two 
beams.  Figure  3  shows  a  schematic  of  a  typical  FTS.  When  the  fixed  and  movable 
mirrors  are  an  equal  distance  from  the  beam  splitter,  the  OPD  is  zero  and  therefore 
constructive  interference  will  occur.  As  the  mirror  moves  and  the  OPD  changes,  the  two 
beams  will  move  in  and  out  of  phase  causing  constructive  and  destructive  interference 
patterns.  In  this  way,  an  interferogram  can  be  constructed  that  contains  the  light  intensity 
as  a  function  of  OPD.  For  monochromatic  light  this  is  simply  a  sinusoid.  If  multiple 
frequencies  are  contained  in  the  incoming  light,  then  an  interferogram  that  is  a 
summation  of  each  individual  interferogram  for  each  individual  wavelength  is  produced. 
The  easiest  way  to  visualize  this  is  by  looking  at  bi-chromatic  light.  In  Figure  4,  two 
discrete  wavelengths  are  incident  upon  the  interferometer  each  producing  its  own 
interferogram.  Adding  the  two  individual  interferograms  together  gives  the  actual 
interferogram  output  from  the  interferometer. 
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Figure  3.  Typical  FTS  layout.  Light  enters  the  system  from  the  scene  on  the  right  and  then  is  split 
by  the  beam  splitter.  Each  path  then  is  reflected  back  to  the  beam  splitter  off  of  the  retro  reflectors 
and  therefore  recombined  before  hitting  the  detector.  As  one  of  the  mirrors  moves  the  phase 
difference  between  the  recombined  light  causes  an  interference  pattern  called  an  interferogram. 
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Figure  4.  The  picture  on  the  left  is  the  frequency  spectrum  of  the  incoming  light  showing  two 
different  wavelengths.  On  the  right  is  the  interferogram  produced  by  this  bi -chromatic  light.  The 
higher  frequency  lower  intensity  wave  is  added  to  the  lower  frequency  higher  intensity  wave  to 
produce  the  resultant  interferogram  which  is  a  superposition  of  the  two  individual  waves. 


The  same  thing  holds  true  for  a  broad-band  spectrum.  Each  spectral  component 
of  the  light  source  produces  its  own  interferogram  with  its  characteristic  period,  and 


whose  amplitude  is  weighted  by  the  relative  spectral  intensity.  An  interferogram  such  as 
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this  is  shown  in  Figure  5.  At  zero  path  difference  (ZPD)  all  wavelengths  constructively 


interfere  which  causes  the  characteristic  center  burst  shown.  At  maximum  path 
difference  (MPD)  the  interference  largely  cancels  out  the  signal. 
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Figure  5.  A  broad-band  interferogram  generated  by  the  Telops  Hypercam. 

One  nice  thing  about  an  FTS  is  that  the  spectral  resolution  is  proportional  to 
MPD.  It  is  common  to  define  the  spectral  resolution  of  an  FTS  as 

1 


8<j  = 


and  the  resolving  power  as 


4  d 


S<7  (7 


(2) 


(3) 


where  d  is  the  total  displacement  of  the  mirror  from  ZPD  [2],  This  allows  the  user  to 

tailor  the  spectral  resolution  to  the  application. 

In  the  ideal  case,  to  produce  a  spectrum  from  the  measured  interferogram  one 

must  compute  the  inverse  cosine  Fourier  transform  of  the  intensity.  The  superposition  of 

waves  in  the  interferogram  with  the  mean  value  subtracted  out  is  given  by 
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(4) 


I(x)  =  B(<j)CO$,(2n<7X)d<J 

where  I(x )  is  the  intensity  of  the  interferogram  as  a  function  of  OPD,  x,  and  B(cr)\s  the 
intensity  of  the  light  as  a  function  of  wavenumber,  <J  .  The  spectrum  can  then  be 
recovered  by  the  inverse  Fourier  transform  given  by 

poo 

B((j)-\  I(x)cos(27T(Tx)dx .  (5) 

J  -00 

In  reality,  the  limits  of  the  integration  in  Equations  (4)  and  (5)  are  determined  by  the 
physical  constraint  on  the  distance  the  mirror  can  travel  and  will  be  some  finite  value. 
Additionally,  the  cosine  transform  is  used  in  the  ideal  case  because  the  symmetry  of  the 
interferograms  results  in  a  real-valued  Fourier  transform.  In  real  world  measurements, 
the  interferograms  are  not  exactly  symmetric  and  some  phase  shifts  result  in  a  complex 
interferogram  and  therefore  a  complex  Fourier  transform  and  a  complex  spectrum.  The 
interferogram  is  also  measured  in  discrete  increments  and  therefore  requires  the  use  of  a 
discrete  Fourier  transform  (DFT)  as  opposed  to  a  continuous  one.  Taking  these  factors 
into  consideration,  Equation  (5)  becomes 

=  (6) 

Nft 

where  N  is  the  number  of  discrete  samples  and  x  .  is  the  distance  of  the  mirror  from 

ZPD  at  those  sample  points.  In  the  case  of  the  Telops,  the  distance  between  sampling 
points  is  determined  by  the  HeNe  laser  shown  in  Figure  3.  Since  the  HeNe  laser  is  a 
monochromatic  light  source,  its  interferogram  is  a  sinusoid  whose  wavelength  is  the 
wavelength  of  the  HeNe  laser  light.  The  interferogram  is  sampled  every  three  HeNe 
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wavelengths.  In  order  to  meet  the  Nyquist  sampling  criteria  for  the  irradiance  spectrum, 
the  irradiance  spectrum  must  be  band  limited  to  a  maximum  wavenumber  of 


cr 


max 


1 

4A 


(7) 


or  a  minimum  wavelength  of 


2,„  =  4A 


(8) 


where  A  is  the  distance  between  sampling  points.  Since  the  HeNe  wavelength  is  632.8 
nm,  that  gives  a  sampling  distance  of  A  =  1.898  xlO-4  cm  which  gives  cr^  =  1317  cm-1 
or  2min  =  7.6  pm.  The  detector  itself  has  good  sensitivity  over  the  8-12  pm  band  [3], 

There  are  advantages  and  disadvantages  to  using  an  FTS.  In  atomic  spectroscopy 
the  line  intensities  often  vary  by  orders  of  magnitude  from  line  to  line.  The  ratio  between 
the  strongest  and  weakest  lines  distinguishable  by  the  instrument  is  called  its  dynamic 
range.  The  FTS  is  typically  superior  at  measuring  both  absolute  and  relative  line 
intensities.  For  an  FTS,  observing  dynamic  ranges  of  1000  is  routine  and  with  care  a 
dynamic  range  of  30,000  can  be  achieved  [4],  As  mentioned  above,  spectral  resolution 
for  an  FTS  is  determined  by  the  maximum  OPD.  This  is  advantageous  because  the  user 
can  control  this  parameter  and  therefore  can  control  the  spectral  resolution.  This  prevents 
unnecessary  oversampling  and  reduction  in  the  signal -to-noise  ratio  (SNR).  The  FTS 
also  benefits  from  the  multiplex  advantage,  that  is,  it  collects  all  wavelengths  at  once  on 
every  pixel  as  opposed  to  spatially  splitting  the  signal  into  spectral  bands.  This  provides 
an  SNR  boost  compared  with  dispersive  spectrometers  although  this  is  offset  by  the  fact 
that  shot  noise  is  also  captured  for  the  entire  spectral  band  as  well.  A  disadvantage  for 


FTS  based  remote  sensing  lies  in  the  fact  that  it  takes  time  to  move  the  mirror  through  its 
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complete  range  of  OPDs.  This  is  problematic  for  observing  scenes  that  are  very  dynamic. 
Additionally,  there  can  be  issues  if  the  FTS  is  mounted  on  a  moving/vibrating  platform 
since  it  is  sensitive  to  any  misalignments. 

Polarimetric  Imaging  in  the  LWIR 

Introduction 

Polarimetric  imaging  measures  the  polarization  state  of  light  coming  from  the 
scene  being  imaged.  This  information  can  help  classify  materials  and  identify  objects  of 
interest  for  remote  sensing  and  military  applications.  At  wavelengths  shorter  than  3  pm , 
the  polarization  signature  of  objects  in  the  scene  is  dominated  by  surface  reflections  and 
scattering.  In  the  LWIR  band,  however,  objects  in  the  scene  are  self-emitting  measurable 
radiation  proportional  to  the  objects  temperature  and  emissivity.  This  self-emitted 
radiation  from  the  material  can  become  partially  polarized  upon  transmission  across  the 
surface  boundary  [5],  In  the  LWIR,  it  is  therefore  common  to  have  contributions  from 
both  the  reflected  and  emitted  radiance  from  an  object.  This  poses  a  difficulty  since  the 
reflected  polarization  component  is  orthogonal  to  the  emitted  polarization  component. 
These  two  orthogonal  polarization  states  add  incoherently  and  result  in  a  reduced 
apparent  degree  of  polarization  [6],  For  this  reason,  the  atmospheric  conditions  play  a 
large  role  in  polarimetric  remote  sensing  applications  in  the  LWIR.  For  a  good 
polarimetric  signature,  there  needs  to  be  thermal  contrast  between  the  surrounding 
atmosphere  which  is  reflecting  off  of  the  target  and  the  target  itself.  If  these  two 
contributions  are  equal,  then  the  apparent  degree  of  polarization  goes  to  zero.  These 
effects  were  demonstrated  by  Felton  et  al.  [7]  who  studied  polarimetric  imagery  during  a 
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multi-day  diurnal  cycle.  They  found  that  in  the  LWIR,  cloud  cover  reduced  the 
polarimetric  contrast  observed  from  their  target,  which  was  the  hull  of  a  tank.  In  the 
LWIR,  clouds  are  similar  to  blackbodies  that  are  warm  in  comparison  to  the  background 
of  a  clear  sky  where  the  cold  upper  atmosphere  is  visible.  They  showed  that  the  contrast 
was  reduced  to  the  noise  level  during  totally  overcast  conditions.  It  is  noted  that  the 
tank’s  engine  was  not  running  so  the  hull  was  allowed  to  drift  with  atmospheric 
temperature.  It  is  also  interesting  to  note  that  under  the  same  conditions,  polarimetric 
contrast  increased  in  the  MWIR,  which  indicates  the  increased  self-emission  component 
in  the  LWIR  versus  the  reflected  component  in  the  MWIR. 

Stokes  Vector  Representation  of  Polarized  Light 

It  is  convenient  to  represent  the  polarization  state  of  light  as  a  Stokes  vector.  In 
1852,  Stokes  showed  that  fully  polarized  light  can  be  represented  as  [8] 

(4c  +  4,  )2  =  (4,  -  4,  )2  +  (2S0xS0 y  COS  4>f  +  (2£0x£0 y  sin  (9) 


where  s0x  and  s0y  are  the  amplitudes  of  the  x  and  y  components  of  the  electric  field 


and  <j)  is  the  phase  difference  between  the  x  and  y  components  of  the  electric  field. 
The  Stokes  parameters  are  then  defined  based  on  this  formulation  where 


*0  £0x  +  £0y 

si  =  4,-4, 

s2  =  2s0x£0y  COS  (f) 

S3  =  2£0x£0y  S^n  ^ 


(10) 


and 


2  2  2  2 

S0  =  S1  +  ^2  +  ^3  ■ 


(11) 
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The  Stokes  parameters  are  real  observable  quantities  that  can  be  measured  by  a  detector 
in  the  form  of  optical  intensities  or  radiometric  energies.  The  s0  term  describes  the  total 

energy  in  the  field.  The  s,  term  describes  the  amount  of  horizontal  or  vertical  linear 
polarization.  The  s2  term  describes  the  amount  of  ±45  degree  linear  polarization  and  the 
^  term  describes  the  amount  of  left  or  right  handed  circular  polarization.  In  terms  of  the 
total  scene  radiance,  Ls ,  these  Stokes  parameters  can  be  arranged  into  a  column  vector 
such  that 


5  = 


Lh  ~  Lv 

L45  ~  A35 

V  LKCp  ~  LLCp  j 


(12) 


whre  Lh  and  Ly  are  the  horizontal  and  vertical  components,  L45  and  Ll35  are  the  45 
degree  and  135  degree  components,  and  LRCP  and  Llcp  are  the  right-handed  and  left- 
handed  circular  polarization  components.  The  Stokes  vector  is  often  normalized  by 
dividing  each  term  by  the  s0  element  so  that 


5  = 


^o=l" 


V  S3  J 


(13) 


where  each  term  now  represents  the  normalized  Stokes  parameter.  In  this  form  s, ,  s2 , 


and  s3  range  from  -1  to  1  with  the  constraint  that  s0  >  ^sf  ~+si  +s32  .  Positive  values 
indicate  horizontal,  +45°,  or  RCP,  whereas  negative  values  indicate  vertical,  135°,  or 
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LCP.  This  gives  a  simple  way  to  represent  the  polarization  state  of  a  beam  as  illustrated 
in  Figure  6. 


Figure  6.  Normalized  Stokes  vectors  used  to  represent  the  polarization  state  of  light. 

Measuring  the  Stokes  parameters  involves  imaging  the  scene  with  the  polarizer 
oriented  at  various  angles.  Since  the  Telops  implements  a  division  of  time  method  for 
collecting  at  different  polarizer  angles,  it  is  important  that  the  lighting  conditions  and 
positioning  of  objects  in  the  scene  and  the  position  of  the  instrument  remain  as  constant 
as  possible  for  a  given  set  of  polarized  images.  A  common  method  of  determining  the 
Stokes  parameters  for  linearly  polarized  light  is  known  as  the  Modified  Pickering 
Method.  With  this  method,  the  incoming  light  is  measured  with  the  polarizer  oriented  at 
0,  45,  90,  and  135  degrees.  The  Stokes  parameters  are  then 

( Lq  +  T45  +  +  Z^35 ) 

2 

s=  •  (14) 

^45  ~  A  35 
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Once  the  Stokes  parameters  are  measured,  they  can  then  be  used  to  calculate  the  degree 
of  linear  polarization  (DoLP)  and  the  angle  of  polarization  (AoP)  as  follows: 


DoLP  = 


A2 +4 


0 


(15) 


]_ 

2 


AoP  =  —  tan  1 1  2 


\si  J 


(16) 


It  should  be  noted  that  in  passive  remote  sensing,  s3  is  usually  close  to  zero  so  it  is  only 

important  to  characterize  the  linear  polarization  state  of  the  beam  [9].  Since  the  polarizer 
on  the  Telops  is  a  linear  polarizer,  this  will  be  the  case  for  this  research.  Some  nice 
features  of  the  Stokes  vectors  are  that  they  work  for  partially  polarized  light  and  they  are 
additive  in  nature  as  long  as  the  radiation  is  incoherent.  The  combined  radiance  and 
polarization  from  multiple  sources  of  arbitrary  polarization  is  simply  the  addition  of  their 
Stokes  vectors.  For  example,  the  combination  of  horizontally  and  vertically  polarized 
light  is  given  by 


S  -  SH+Sy~ 


f 1 ) 

f  n 

f l ) 

1 

4  + 

-i 

4  = 

0 

2  L 


so  that  the  result  is  unpolarized  light.  The  combination  of  horizontal  and  45  degree 
polarized  light  is 


S  -SH+S 45  - 


(l) 

(l) 

f  l  ) 

l 

4  + 

0 

4  = 

0.5 

w 

[0.5) 

2L  . 
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Mueller  Matrix  Representation 

The  interaction  of  Stokes  vectors  with  their  environment  is  modeled  with  Mueller 
matrices.  A  Mueller  matrix  converts  an  input  Stokes  vector  into  an  output  Stokes  vector, 
i.e. 


S out  -  M  •  Sin  - 


moo 

m01 

m02 

m03^ 

"h 0 

"hi 

"hi 

mi3 

m20 

m2j 

m22 

m23 

"ho 

m31 

m32 

m33  J 

r  o  \ 


\SiJ 


(17) 


Any  optical  element  can  be  represented  by  its  Mueller  matrix.  For  light  traveling  through 
an  optical  system  consisting  of  many  optical  elements,  a  combined  system  Mueller 
matrix  can  be  found  through  the  matrix  multiplication  of  each  of  the  individual  elements’ 
Mueller  matrices.  In  this  case  the  output  Stokes  vector  is  given  by 

(18) 

where  is  the  first  optical  element  with  which  the  incoming  Stokes  vector  interacts 


and  M„  is  the  last.  The  Mueller  matrix  of  an  ideal  polarizer  oriented  at  an  angle,  6p , 


measured  positive  when  rotating  counterclockwise  from  the  x  -axis  looking  into  the  beam 


is  given  by 


M«9,)  =  - 


1 

cos  26 p 
sin  26 p 
0 


cos  26 p 
cos2  26 p 
sin  26 p  cos  26 p 
0 


sin  26 p  (T 

sin  26 p  cos  26 p  0 

sin2  26 p  0 

0  0y 


(19) 


No  polarizer  is  perfectly  ideal  and  there  are  two  common  metrics  that  characterize  how 
close  to  ideal  a  linear  polarizer  is.  These  are  the  extinction  ratio  given  by 


19 


ER  =  (20) 

T  ■ 

min 

where  rmax  and  rmin  are  the  maximum  and  minimum  transmission  of  linearly  polarized 
light  through  the  polarizer,  and  the  diattenuation  which  is  given  by 


D  = 


T  —  T 

max _ min 

T  4-7" 

max  min 


(21) 


Diattenuation  has  the  benefit  of  ranging  from  1  for  an  ideal  polarizer  to  0  for  a  non¬ 
polarizing  element.  Taking  diattenuation  into  consideration,  the  Mueller  matrix  for  a 


linear  polarizer  becomes 


\_ 

2 


q  +  r 

(q  -  r)  cos  2  6 
(q-r)  sin  26  p 


V 


0 


(q  -r)cos26p 

(q  +  r)  cos 2  26 p  +  2 <Jqr  sin 2  2 6 p 
(q  +  r- 2 -Jqr) sin  20 p  cos 26 p 
0 


(q  -  r)  sin  26 p 

(q  +  r- 2 s[qr) sin  26 p  cos  26 p 
(q  +  r) sin2  2 6 p  +  2^qr  cos2  2 6 p 
0 


0 

0 


2  Jqr) 


(22) 


where  q  =  rmax  and  r  =  rmjn .  The  Stokes  vector  and  Mueller  matrix  approach  will  be 


used  heavily  in  this  thesis  when  developing  the  polarimetric  system  model  and  calibration 
methodology. 


Instrumentation 

The  instrument  used  for  this  research,  shown  in  Figure  7,  is  an  LWIR  Field 
Portable  Imaging  Spectrometer  (FIRST)  developed  by  Telops,  Inc.  known  as  a  Hyper¬ 
cam.  It  has  been  modified  with  a  rotatable  zinc  selenide  (ZnSe)  linear  wire-grid  polarizer 
on  the  front  end  to  add  the  polarimetric  capability.  The  spectrometer  inside  the 
instrument  is  essentially  a  Michelson  interferometer  that  uses  corner-cubes  instead  of  flat 
mirrors  to  reduce  the  effect  of  mechanical  vibrations  as  the  mirror  is  scanning.  The 
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imaging  capability  is  made  possible  by  a  320  x  256  photovoltaic  focal  plane  array  (FPA) 
made  of  Mercury  Cadmium  Telluride  (HgCdTe)  with  a  detector  pitch  of  30  //m .  The 
spectral  resolution  can  be  tuned  by  adjusting  the  distance  the  scan  mirror  travels.  This 
provides  a  spectral  resolution  ranging  from  0.25  cm"1  to  150  cm"1.  The  field  of  view 
(FOV)  of  the  instrument  can  also  be  controlled  electronically  by  windowing  down  to  any 
sub-region  of  the  array.  The  time  it  takes  to  collect  a  full  interferogram  is  proportional  to 
the  spectral  resolution  and  the  frame  rate  is  proportional  to  the  FOV,  so  the  ability  to 
adjust  both  of  these  parameters  allows  more  flexibility  in  achieving  collection  goals.  The 
maximum  frame  rate  at  full  window  is  300  fps  and  for  a  128  x  128  sub-region  it  increases 
to  1400  fps.  At  1400  fps  the  instrument  can  collect  a  full  data  cube  at  4  cm"1  resolution 
every  two  seconds.  The  ZnSe  wire  grid  polarizer  is  housed  in  a  rotation  stage  that  allows 
±180°  rotation  and  achieves  an  orientation  accuracy  of  less  than  0.1°  and  an  orientation 
stability  of  less  than  0.01°.  Additional  parameters  are  provided  in  Table  1. 
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Figure  7.  The  Telops  P-HSI  system.  The  picture  on  the  bottom-middle  shows  the  interferometer 
module  and  the  picture  on  the  bottom-right  shows  the  polarization  module.  (Image  courtesy  of 

Telops,  Inc.). 

Table  1.  Telops  P-HSI  system  parameters 


Parameter 

Specification 

Effective  focal  length 

8.6  cm 

Aperture  diameter 

4.3  cm 

F /# 

2.0 

Detector  array  size 

320  x  256 

Detector  pitch 

30  jum 

IFOV 

0.35  mrad 

FOV 

5.1  x  6.4  degrees 

Spectral  range 

8-12  jum 

Spectral  resolution 

0.25-150  cm'1 

Number  of  spectral  channels 

3  to  1600 

Maximum  frame  rate 

300  Hz  (320  x  256) 

1400  Hz  (128  x  128) 

Polarizer  type 

Wire-grid  (linear) 

Polarizer  material 

ZnSe 
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Polarizer  orientation  accuracy  <0.1° 

Polarizer  orientation  stability  <  0.01° 


Previous  Work  Related  to  Spectro-Polarimetric  Performance  and  Calibration 

This  section  will  introduce  several  previous  works  that  are  relevant  to  this  thesis. 
As  described  above,  the  Telops  instrument  is  composed  of  many  different  technologies 
which  could  each  have  their  own  section.  This  section,  however,  will  focus  on  the 
polarimetric  aspect  and  the  calibration  aspect  in  general. 

This  section  would  not  be  complete  without  the  mention  of  Sir  George  Stokes  and 
Hans  Mueller.  In  1852,  Stokes  published  a  paper  called  “On  the  composition  and 
resolution  of  streams  of  polarized  light  from  different  sources”  [8]  which  was  published 
by  the  Cambridge  Philosophical  Society.  This  paper  was  foundational  to  the  way 
polarized  light  is  represented  today  with  so  called  Stokes  vectors.  How  to  describe  the 
interaction  of  this  polarized  light,  in  Stokes  vector  representation,  with  the  surrounding 
world  was  developed  by  Hans  Mueller  in  1943.  In  [10],  he  describes  a  4  x  4  matrix 
which  converts  an  input  Stokes  vector  to  an  output  Stokes  vector  after  having  been 
reflected  or  transmitted  through  a  material.  The  principles  pioneered  by  these  two 
scientists  are  used  heavily  in  this  thesis.  Much  of  the  theory  and  applications  of  polarized 
light,  Stokes  vectors,  and  the  Mueller  calculus  was  later  covered  in  the  texts  written  by 
Collet,  Goldstein,  Chipman,  and  Hecht  [11,  12,  13,  14,  15].  An  important  aspect  of 
polarimetry  in  the  LWIR  is  the  concept  of  emission  polarization.  The  background 
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required  to  understand  the  phenomenon  of  emission  polarization  is  provided  by  Sandus  in 
his  1965  journal  article  in  Applied  Optics  [5]. 

The  use  of  polarized  light  in  remote  sensing  and  imagery  applications  started  to 
gain  ground  as  early  as  the  1980s  with  works  such  as  “Polarization  Imagery”  by 
Walgraven  [16]  and  “Polarization  Imaging”  by  Solomon  [17]  both  published  in  1981. 

The  field  has  continued  to  grow  ever  since.  More  recent  work  used  by  the  author  to 
develop  an  understanding  of  polarimetry  in  remote  sensing  applications  was  Schott’s 
books,  “Fundamentals  of  Polarimetric  Remote  Sensing”  and  “Remote  Sensing,  the  Image 
Chain  Approach”  [9,  18].  With  regards  to  hyperspectral  remote  sensing,  Eismann’s  book 
titled  “Hyperspectral  Remote  Sensing”  was  very  useful;  in  particular  the  chapters, 
“Fourier  Transform  Spectrometer  Design  and  Analysis”  [2]  and  “Spectrometer 
Calibration”  [19]  were  quite  useful  for  this  work. 

For  the  specific  application  of  polarimetric  calibration,  there  are  many  useful 

sources  in  which  to  turn.  The  work  that  is  the  basis  of  the  experimental  polarimetric 

calibration  effort  completed  in  this  thesis  was  written  by  Persons  et  al.  and  is  titled  “A 

proposed  standard  method  for  polarimetric  calibration  and  verification”  [20].  In  their 

work,  they  describe  a  method  for  using  linear  combinations  of  Stokes  vectors  to  generate 

highly  accurate  input  states.  From  these  input  states  and  the  measurement  of  the 

calibrated  instrument  response,  the  data  reduction  matrix  can  be  generated.  Furthermore, 

they  describe  how  to  present  the  data  in  such  a  way  that  is  useful  in  understanding  the 

non-ideal  effects  of  the  system.  According  to  them,  the  method  described,  which  allows 

direct  measurement  of  the  system  effects,  is  an  improvement  over  other  methods  which 

rely  on  comparing  modeled  calibration  parameters  with  measured  data  since  it  is  difficult 
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to  model  higher  order  effects.  Kudenov,  Pezzaniti,  and  Gerhart  utilize  similar  methods 
and  provide  further  mathematical  detail  and  modeling  as  well  as  experimental  validation 
in  their  2009  paper  on  their  microbolometer-infrared  imaging  Stokes  polarimeter  [21]. 

The  papers  written  by  Smith  et  al.  [22]  in  2000  and  Blumer  et  al.  [23]  in  2002 
both  discuss  the  calibration  of  the  Multispectral  Infrared  Stokes  Imaging  Polarimeter. 

The  first  paper  discusses  calibration  in  the  MWIR  and  the  later  paper  discusses 
calibration  in  the  LWIR.  In  both  cases,  the  authors  found  that  there  was  a  variation  in 
instrument  response  as  a  function  of  retarder  orientation.  In  the  first  paper,  the  problem 
was  reduced  by  an  order  of  magnitude  by  tilting  the  polarizer  angle  slightly  to  reduce  the 
narcissus  effect.  The  effect,  however,  was  never  fully  calibrated  out.  In  the  second  paper, 
a  different  method  that  involves  performing  a  non-uniformity  correction  (NUC)  at  each 
retarder  angle  was  used  to  achieve  much  better  results.  The  current  calibration  method 
for  the  Telops  also  involves  calibration  at  each  polarizer  angle.  It  will  be  shown  in  the 
Polarimetric  Calibration  chapter  that  this  removes  any  polarimetric  biases  of  the  system 
to  the  extent  that  the  instrument  polarizer  is  ideal. 

A  more  recent  development  for  LWIR  imaging  polarimeters  includes  the  micro¬ 
grid  polarizer  array.  This  involves  dividing  the  FPA  into  a  grid  of  super-pixels  where 
each  super-pixel  consists  of  four  pixels,  each  with  a  polarizing  element  placed  in  front  of 
it  at  0,  45°,  90°,  and  135°  respectively.  By  using  a  division  of  focal  plane  (DoFP)  scheme 
rather  than  a  division  of  time  instrument  such  as  the  Telops,  the  user  gets  near  real-time 
polarization  imagery  with  no  thermal  drift  in-between  measurements.  Any  errors  in  the 
NUC,  however,  will  cause  false  polarization  signatures  similar  to  thermal  drift  in  a 

division  of  time  instrument.  The  papers  by  Bowers  et  al  in  2006  [24]  and  2008  [25] 
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discuss  the  NUC  of  the  instrument  and  provide  useful  insight  into  the  effect  these  have  on 
the  polarimetric  results,  while  Hubbs  et  al.  [26]  discuss  the  radiometric  and  polarization 
characteristics  of  the  micro-grid  FPA  in  their  2006  paper. 

A  couple  of  papers  provided  useful  information  with  regards  to  determining  and 
understanding  the  noise-equivalent  spectral  radiance  (NESR)  and  the  noise-equivalent 
spectral  degree  of  linear  polarization  (NESDoLP).  The  first  is  a  paper  written  by  Telops 
Inc.  in  2006,  which  describes  a  model  to  predict  the  NESR  of  the  Telops  instrument  and 
compares  this  to  experimental  results  [3].  The  experimentally  determined  NESR  in  the 
paper  is  consistent  with  the  results  found  in  this  thesis.  The  second  paper,  written  by 
Jones  and  Persons  [27]  in  2007,  develops  a  model  to  predict  the  performance  of  micro¬ 
grid  polarimeters.  The  polarimetric  performance  is  described  in  terms  of  the  noise- 
equivalent  DoLP  (NEDoLP).  The  results  of  the  model  provide  insight  into  how  various 
system  parameters  affect  the  NEDoLP,  and  hence,  provide  a  better  understanding  of  the 
results  obtained  in  this  thesis. 

Many  of  the  papers  referenced  in  this  thesis  and  the  techniques  that  are  derived 
from  them  are  done  on  a  broadband  imaging  basis  and  not  a  hyperspectral  basis.  Despite 
this,  it  was  not  too  difficult  to  extrapolate  the  techniques  obtained  from  these  papers  to 
the  hyperspectral  case  by  treating  each  waveband  as  an  individual  image.  Thus,  for 
example,  rather  than  obtaining  a  2-D  array  of  data  reduction  matrices  (one  for  each 
pixel),  a  3-D  cube  was  generated  with  the  data  reduction  matrix  for  each  pixel  at  each 
waveband. 
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III.  Standard  HSI  Calibration 


This  chapter  discusses  the  calibration  procedures  which  must  be  accomplished  for 
the  Telops  HSI  instrument  whether  the  polarization  module  is  present  or  not.  Both 
sections  discuss  the  theory  that  supports  the  calibration  procedures  and  then  provide 
experimental  data  demonstrating  the  effectiveness  of  the  calibrations.  The  insight 
developed  in  the  first  section,  which  discusses  FPA  non-uniformity  and  radiometric 
calibration,  will  be  particularly  useful  later  on  when  discussing  the  polarimetric 
calibration. 

FPA  Non-uniformity  and  Radiometric  Calibration 

When  collecting  radiometric  data,  the  measured  signal  at  the  detector  includes  a 

contribution  from  both  the  scene  as  well  as  an  offset  radiance  generated  from  the 

instrument  itself.  This  is  particularly  evident  in  the  thermal  IR  where  optical  components 

in  the  system  as  well  as  electronics  are  all  thermally  emitting  photons.  This  signal, 

consisting  of  scene  photons  and  instrument  photons,  is  then  subjected  to  transmission 

losses  in  the  optical  system  as  a  function  of  wavelength  as  well  as  the  spectral  response  of 

the  detector  elements.  The  photons  must  then  generate  electrons  to  be  stored  by  the  pixel 

elements,  read  out  by  the  detector  electronics  and  then  ultimately  converted  to  digital 

numbers.  These  digital  numbers  must  then  be  converted  to  useful  radiometric  units 

through  a  radiometric  calibration  process.  Due  to  pixel  nonuniformity,  this  radiometric 

calibration  must  be  accomplished  on  a  pixel-by-pixel  basis.  If  this  is  not  done,  fixed 

pattern  noise  from  array  nonuniformity  will  dominate  the  signal  rendering  it  useless  as 

will  be  shown  below.  To  remove  this  nonuniformity,  a  two  point  calibration  is 
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performed.  It  should  be  noted  that  the  two  point  calibration  assumes  that  the  focal  plane 
array  (FPA)  exhibits  a  linear  response.  Any  level  of  FPA  response  nonlinearity  will 
result  in  residual  spatial  nonuniformity  even  after  the  two-point  radiometric  calibration 
[19].  In  general,  however,  good  results  have  been  achieved  with  the  Telops  using  a  two 
point  calibration. 

The  two  point  calibration  utilizes  a  spatially  uniform  calibration  source  of  know 
radiance  that  is  first  set  to  a  low  temperature  and  then  to  a  high  temperature  during 
measurement.  The  two  temperatures  are  set  so  that  they  bracket  the  expected  radiance 
that  will  exist  in  the  scene.  It  is  important  that  the  calibration  source  overfills  the  field  of 
view  (FOV)  of  the  instrument  during  calibration  so  that  spatially  uniform  light  is  incident 
upon  all  pixels  of  the  array.  The  linear  equation  used  to  relate  digital  counts  to  scene 
radiance  for  a  particular  pixel  is  given  by 

YJA)  =  G(A)Ls(A)  +  0(A),  (23) 

where  Ym(/ l)  is  the  measured  signal  in  digital  numbers  as  a  function  of  wavelength,  G(A) 
is  the  instrument  gain  as  a  function  of  wavelength,  Ls  (A)  is  the  spectral  radiance  of  the 
calibration  source,  and  O(A)  is  the  instrument  offset  as  a  function  of  wavelength  in 
digital  numbers.  Using  calibration  sources  at  two  different  temperatures  gives  a  system 
of  two  equations  with  two  unknowns  which  can  then  be  solved  for.  The  Telops  has  two 
on-board  black  body  calibration  sources  that  can  be  set  to  different  temperatures  to 
streamline  the  data  collection  process.  For  a  standard  camera,  this  is  a  straightforward 
process,  but  due  to  the  nature  of  the  FTS,  there  are  some  additional  steps  that  must  take 
place  before  the  gain  and  offset  are  solved  for.  This  has  to  do  with  the  fact  that  the  FTS 
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produces  a  discrete  complex  raw  interferogram  at  each  pixel  rather  than  a  real  raw 


spectrum.  The  imaginary  part  is  a  result  of  phase  errors  introduced  by  asymmetries  in  the 
interferogram  [4]  and  the  discrete  part  is  a  result  of  sampling.  The  discrete  complex  raw 
interferograms  must  then  be  transformed  to  spectra  using  the  complex  discrete  Fourier 
transform  (DFT)  given  by 

YW  =  •)/(*,.)/ ^  ,  (24) 

N 


where  J(xj )  is  the  raw  interferogram,  x  is  the  optical  path  difference,  and  W(Xj)  is  an 

optional  apodization  function  to  reduce  ringing  in  the  spectrum.  The  complex  offset  and 
gain  can  then  be  solved  for  using  the  complex  version  of  equation  (23)  for  the  two 
different  black  body  measurements.  This  gives: 


G(1)  =  M), 
4U)-AU) 


(25) 


and 


(26) 


The  imaginary  part  of  the  raw  spectrum,  gain,  and  offset  is  an  indication  of  noise 
and  can  be  used  as  a  quick  upfront  diagnostic  of  the  quality  of  the  data.  For  example, 
Figure  8  shows  the  real  and  imaginary  parts  for  the  data  set  that  will  be  discussed  below. 
Notice  that  the  imaginary  part  of  each  curve  is  close  to  zero  and  the  real  part  is  much 
higher.  This  is  an  early  indication  that  there  is  a  good  SNR. 
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Figure  8.  The  complex  raw  spectrum  (top),  complex  gain  (middle),  and  complex  offset  (bottom)  for 

the  data  collection. 


A  series  of  images  will  now  be  presented  that  walk  through  the  process  described 
above  to  go  from  a  raw  interferogram  data  cube  to  a  calibrated  spectral  data  cube.  Figure 
9  shows  the  scene  that  was  captured  with  the  Telops  that  will  be  used  for  this  section. 

The  image  was  taken  by  a  visible  camera  mounted  to  the  top  of  the  Telops  that  is  bore 
sighted  with  the  infrared  camera  in  the  FTIR.  The  blue  box  roughly  outlines  the  scene, 
although  it  needs  to  be  shifted  up  and  to  the  right  several  pixels  to  represent  the  true 
scene. 
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Figure  9.  Visible  imagery  taken  of  the  scene  with  a  camera  mounted  on  top  of  the  Telops.  The  blue 
box  roughly  outlines  the  area  of  the  scene  collected  by  the  infrared  camera  in  the  FTIR. 

Figure  10  shows  the  raw  band-averaged  image  as  well  as  the  spectral  information 
contained  within  an  arbitrary  pixel  before  and  after  they  are  converted  from  raw 
interferograms  to  raw  spectra.  The  figure  demonstrates  the  importance  of  the 
nonuniformity  correction.  The  fixed  pattern  noise  of  the  array,  before  correction, 
overwhelms  the  signal  and  washes  out  the  image  as  well  as  any  resolvable  spectral 
features.  To  have  a  useful  image,  the  two  point  calibration  discussed  above  must  be 
implemented  to  remove  gain  and  offset. 
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Raw  Spectral  Image 


Raw  Interferogram  Image 
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Figure  10.  On  the  right,  a  raw  band-averaged  interferogram  image  is  shown  along  with  an 
interferogram  from  a  random  pixel.  On  the  left  a  band-averaged  image  is  shown  after  having  been 
converted  from  raw  interferograms  to  raw  spectra.  The  raw  spectrum  from  a  pixel  is  shown  on  the 

bottom. 


Figure  1 1  shows  the  result  of  applying  the  radiometric  calibration.  The  spectrum 
for  the  same  pixel  is  displayed  next  to  the  band-averaged  image.  For  this  data  set  the 
resolution  was  approximately  8  cm'1.  At  this  resolution,  several  water  absorption  features 
as  well  as  the  ozone  absorption  feature  at  1043  cm'1  are  now  coarsely  resolved  revealing 
that  this  pixel  is  one  of  the  sky  pixels  above  the  building.  While  this  is  a  calibrated 
image,  it  is  still  unprocessed  and  many  bad  pixels  can  be  seen.  The  image  is  also  flipped 
about  the  vertical  axis.  Figure  12  shows  the  fully  calibrated  and  processed  band-averaged 
image.  Here,  the  bad  pixels  have  been  removed,  the  image  has  been  oriented  properly, 
and  a  different  color  map  is  used.  With  an  image  that  contains  a  high  dynamic  range 
such  as  this  one,  any  residual  FPN  is  negligible  and  can  no  longer  be  distinguished  in  the 
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image.  The  radiometric  calibration  principles  discussed  in  this  section  will  be  useful  in 


understanding  the  full  polarimetric  calibration  which  will  be  discussed  later. 

Calibrated  Unprocessed  Radiance  Image 


— 

Figure  11.  Calibrated  band-averaged  radiance  image  (left)  and  a  spectrum  (right)  for  one  pixel 
looking  at  the  sky  above  the  building.  Several  water  absorption  features  are  coarsely  resolved  as  well 

as  the  ozone  feature  at  1043  cm'1. 
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Figure  12.  Calibrated  and  processed  band-averaged  image  of  the  scene.  The  bad  pixels  have  been 
removed,  the  image  has  been  oriented  properly,  and  a  more  natural  color  map  is  used. 
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Off-Axis  Spectral  Calibration 

Spectral  calibration  of  the  Telops  FTS  is  made  possible  with  the  use  of  a  HeNe 
laser.  The  HeNe  laser  provides  a  monochromatic  source  at  632.816  nm  that  results  in  a 
sinusoidal  interferogram.  If  the  zero  optical  path  difference  location  is  known,  the  phase 
of  the  HeNe  interferogram  can  be  used  to  precisely  measure  the  optical  path  difference  as 
the  mirror  position  changes,  which  in  turn  allows  for  the  precise  measurement  of  spectral 
features  when  performing  the  FFT.  This  method  only  works  precisely,  however,  when 
the  light  from  the  scene  travels  the  exact  path  through  the  instrument  as  the  HeNe  laser 
light  travels.  For  a  pixel  that  is  off-axis  with  respect  to  the  HeNe  laser,  the  optical  path 
difference  traveled  by  light  reaching  that  pixel  is  slightly  different  than  that  of  the  HeNe 
laser.  Because  of  this,  the  spectral  axis  for  pixels  that  are  off-axis  from  the  HeNe  laser  is 
scaled  with  respect  to  the  true  values.  This  results  in  the  necessity  for  an  off-axis  spectral 
calibration  when  highly  accurate  spectral  feature  measurement  is  desired  for  all  pixels. 
The  effect  is  small  since  the  off-axis  angles  spanned  by  the  detector  are  small,  but  for 
highly  accurate  spectral  measurements,  it  still  needs  to  be  corrected. 

An  additional  effect  to  be  considered  for  off-axis  pixels  is  self  apodization  caused 

by  the  spatial  fringe  pattern  created  when  light  travels  through  an  FTS.  This  fringe 

pattern  is  a  function  of  mirror  displacement.  The  further  the  mirror  travels  away  from 

ZPD,  the  higher  the  spatial  frequency  of  the  fringe  pattern.  Since  the  pixels  integrate 

spatially  over  the  fringe  pattern,  if  the  fringe  pattern  causes  a  significant  phase  variation 

over  the  pixel,  it  will  decrease  the  modulation  depth  of  the  interferogram.  The  effect 

increases  with  increasing  mirror  displacement.  The  first  part  of  this  section  provides  an 
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analysis  of  the  apodization  effect  and  the  second  part  of  this  section  discusses  the  off-axis 
spectral  calibration. 

In  an  imaging  FTS,  the  irradiance  on  the  FPA  as  a  function  of  OPD  can  be  written 


as 


poo  1 

E(x,  y,d)  =  JQ  -E0  (x,  y,  a) 


1  +  cos 


^  4/rcrd  ^ 
cos#(x,  y) 


da 


(27) 


where  x  and  y  are  the  distance  of  a  pixel  on  the  FPA  away  from  optical  path  center  and 
d  is  the  OPD  [2],  It  is  the  dependence  of  Equation  (27)  on  9(x,  y)  that  generates  the 
fringe  pattern  of  a  single  frame  when  viewing  a  uniform  input  radiance  over  the  FOV. 
9(x,  y)  is  given  by 


9(x,  y)  =  tan  1 


Jl+7 

f 


(28) 


where  /  is  the  focal  length.  To  analyze  the  off-axis  effects,  it  is  necessary  to  determine 
the  phase  variation  of  the  fringe  pattern  over  the  size  of  one  pixel.  From  Equation  (27), 
the  phase  term  is  given  by 


</>  = 


4  7iad 
cos  9(x,  y)  ’ 


(29) 


taking  the  partial  derivative  of  this  with  respect  to  6  gives  the  differential  phase  change 
with  field  angle,  or 


d(/>  = 


Anad  sin  9 
cos2  9 


d9. 


(30) 


From  Equation  (28)  the  angular  variation  over  the  width  of  a  detector  element,  w ,  is 
given  by 
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(31) 


where  R  is  the  resolving  power  given  by 

R  =  4da.  (33) 

As  mentioned  before,  if  this  phase  variation  over  a  detector  element  is 
significant,  it  causes  apodization  and  thus  decreased  spectral  resolution.  A  phase 
variation  of  2 n  is  similar  to  a  cosine  apodization  where  the  zeroes  of  the  cosine  are 
aligned  with  the  edges  of  the  pixel.  If  In  is  plugged  into  Equation  (32)  for  the  phase 
variation  and  the  field  angle  is  solved  for,  the  result  is 


#  =  sin  1 


(34) 


Therefore,  in  order  to  restrict  the  effect  of  the  phase  variation,  the  field  angle  must  be 
kept  smaller  than  the  value  given  by  Equation  (34).  For  0.25  cm"1  resolution,  which  is 
the  Telops  highest  resolution  setting,  at  a  maximum  wavenumber  of  1250  cm"1, 

R  =  5000 .  This  means  that  the  maximum  field  angle  according  to  Equation  (34)  is 
actually  larger  than  90°  which  is  unphysical.  So  even  if  sin  9  in  Equation  (32)  is  its 
maximum  value  of  1,  a  single  pixel  would  still  not  see  a  2 n  phase  shift.  To  compute  the 
phase  shift  that  a  pixel  will  see,  first  the  maximum  field  angle  must  be  computed  based 
on  the  size  of  the  FPA.  The  Telops  has  a  256x320  array  of  pixels  with  a  pixel  size  of  30 
pm .  The  focal  length  is  86  mm.  This  gives  a  maximum  angle  of 
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Plugging  the  numbers  into  Equation  (32)  gives 


\d(p\  = 


27V  •  5000  •  30//m  •  sin(8. 14° ) 
2  •  86  mm 


=  0.776  rad  =  44.45° 


for  a  single  pixel.  This  will  cause  some  apodization  but  will  still  allow  the  desired 
resolution  to  be  achieved. 


Now  that  the  effect  of  off-axis  apodization  from  the  fringe  pattern  has  been 
quantified  for  the  Telops,  the  off-axis  spectral  scaling  must  be  addressed.  Figure  13 
shows  an  image  that  was  taken  of  the  same  scene  as  in  the  previous  section,  except  at  0.5 
cm'1  spectral  resolution.  The  figure  also  shows  the  results  of  a  type  of  end-member 
analysis  that  utilizes  the  k-means  algorithm.  The  parameter  used  for  distance  in  the  k- 
means  algorithm  was  correlation  and  the  number  of  end-members  was  set  to  10.  The  top 
right  image  groups  end-members  by  color.  What  stands  out  in  this  image  is  the 
concentric  rings  of  color  on  the  front  of  the  uniformly  illuminated  building  and  in  the  sky. 
These  areas  should  have  uniform  intensity  and  spectral  features,  and  yet  the  algorithm 
splits  them  into  different  groups  of  concentric  circles  with  the  center  of  the  circles 
seemingly  in  the  center  of  the  FPA.  This  is  the  result  of  off-axis  spectral  scaling  since 
this  was  not  corrected  in  post-processing.  As  mentioned  before,  the  off-axis  scaling  is  a 
small  effect  but  since  the  algorithm  was  forced  to  find  ten  end  members  in  a  scene  that 
only  has  four  true  spectral  end  members,  it  picked  up  on  the  small  differences  caused  by 
the  scaling.  The  scaling  can  be  confirmed  upon  inspection  of  the  bottom  plot.  The  sharp 
atmospheric  absorption  features  do  not  quite  line  up  for  each  end-member  spectrum. 
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Brightness  Temperature  Map  [K] 
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Figure  13.  Imagery  of  the  same  scene  with  a  0.5  cm'1  spectra  resolution.  Shows  the  results  of  an  end 
member  analysis  that  utilizes  the  k-means  algorithm  with  correlation  as  the  distance  parameter.  The 
top  right  image  groups  the  end-members  by  color  and  the  bottom  plot  shows  the  spectrum  of  each 
end-member  on  a  brightness  temperature  scale.  The  concentric  rings  are  a  result  of  off-axis  spectral 
scaling  which  has  not  been  compensated  for  in  post  processing. 


This  off-axis  scaling  effect  is  typically  compensated  for  by  using  a  modified 
version  of  the  DFT  in  Equation  (24)  for  each  pixel  that  includes  the  6  dependent  phase 
term  in  Equation  (27).  The  corrected  DFT  is 


4  nerd ; 


Y (x,  y,  cr)  =  —  ^ I(x,  y,  d  )e  La,s0<x-y> 

NP 


(35) 


To  implement  this  correction,  optical  path  center  must  be  known  as  precisely  as 
possible.  For  most  laboratory  instruments,  the  optical  path  center  is  known  and  remains 
stable.  For  a  field  portable  instrument  such  as  the  Telops  which  is  often  moved  around 
and  shipped  to  the  field,  the  optical  path  center  is  not  assumed  to  be  constant  and 
therefore  must  be  determined  on  a  periodic  basis.  This  is  accomplished  through  an 
iterative  process  that  utilizes  a  uniform  spectral  calibration  source  with  known  spectral 
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features  that  fills  the  FOV  of  the  detector.  For  the  Telops  LWIR  instrument,  the 
calibration  source  is  a  polystyrene  film  placed  between  the  Telops  and  the  blackbody. 
The  calibration  setup  is  shown  in  Figure  14.  The  polystyrene  film  overfills  the  FOV  of 
the  instrument.  The  polystyrene  film  has  known  spectral  features  within  the  sensitive 
band  of  the  instrument. 


Blackbody 

Polystyrene  film  I 


Figure  14.  Setup  for  generating  known  spectral  features  in  the  LWIR  for  use  during  off-axis  spectral 

calibration. 

The  calibration  scene  must  generally  be  acquired  multiple  times  to  allow  for  data 
cube  averaging  in  order  to  improve  the  SNR  since  the  absorption  lines  must  be  resolved 
on  a  pixel-by-pixel  basis.  Figure  15  shows  the  spectrum  acquired  for  several  pixels  at 
varying  locations  on  the  FPA.  The  box  to  the  right  shows  a  zoomed-in  view  of  the 
absorption  features.  Notice  that  the  feature  does  not  perfectly  line  up  for  each  pixel. 
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Figure  15.  Polystyrene  spectra  of  several  pixels  at  varying  locations  on  the  FPA  before  the  off-axis 
spectral  correction  is  applied.  The  zoomed  in  box  shows  that  the  spectral  features  do  not  perfectly 

line  up  for  each  pixel. 


To  determine  optical  path  center,  an  initial  guess  is  chosen  for  optical  center  and 
random  pixels  are  chosen  from  the  calibration  scene.  The  spectrum  in  these  pixels  is  then 
oversampled  and  divided  by  their  area  for  normalization.  The  pixel  corresponding  to 
optical  center  and  the  instantaneous  field  of  view  (IFOV)  is  then  found  by  minimizing  the 
pair-wise  residual  standard  deviation  between  spectral  features.  Once  optical  center  is 
found,  a  matrix  of  multiplicative  spectral-axis  scale  factors  is  computed  by  determining 
the  1  /  cos  6  term  in  Equation  (35).  The  field  angle  in  Equation  (28)  can  be  re-written  in 
terms  of  pixel  location  on  the  FPA  and  the  instantaneous  field  of  view  (IFOV)  of  the 
pixel  so  that 

0  =  tan"1  ^tan  [(/?  -  )  •  IFOVy  J  +  tan  [(C  -  C0 )  •  IFOVv  f  (36) 

where  R0  is  the  row  index  of  optical  path  center,  C0  is  the  column  index  of  optical  path 

center,  R  is  the  row  index  of  a  given  pixel,  C  is  the  column  index  of  the  pixel,  IFOVy  is 
the  IFOV  in  the  vertical  dimension,  and  IFOVx  is  the  IFOV  in  the  horizontal  dimension. 
The  1  /  cos  0  scale  factor  for  each  pixel  can  then  be  determined  and  placed  in  a  matrix 
given  by 
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(37) 


M(/?,  C)  =  cos 


tan  [(R-Rq)-  IFOVy  J  +  tan  [(C  -  C0 )  ■  1FOV  ]2 


The  matrix  of  scale  factors  can  then  be  used  to  rescale  the  spectral  axis  of  each  off-axis 
pixel  for  any  hyperspectral  data  cube,  thus  allowing  for  more  accurate  spectral 
classification  and  target  detection  of  objects  in  the  scene.  Figure  16  shows  the 
polystyrene  data  after  having  been  corrected  with  the  rescaling  matrix.  Now  the  spectral 
features  at  each  pixel  line  up  almost  perfectly  on  the  spectral  axis. 


Figure  16.  Polystyrene  spectra  of  several  pixels  at  varying  locations  on  the  FPA  after  the  off-axis 
spectral  correction  is  applied.  The  zoomed  in  box  shows  that  the  spectral  features  now  align  almost 

perfectly  along  the  spectral  axis. 


Now  that  the  rescaling  matrix  has  been  determined,  the  end-member  analysis  is 

performed  on  the  corrected  data.  Figure  17  shows  the  results  of  the  same  end-member 

analysis  after  the  off-axis  spectral  scaling  has  been  corrected.  Now  the  absorption 

features  in  the  bottom  plot  all  line  up  properly  and  the  concentric  rings  on  the  front  of  the 

building  and  in  the  sky  due  to  off-axis  scaling  are  gone.  Unlike  before,  the  sky  is  divided 

into  groups  based  on  altitude,  which  is  expected.  What  can  also  be  seen  now  is  the  fringe 

pattern  caused  by  the  interferometer  itself.  This  was  discussed  earlier  and  is 

characterized  by  alternating  lighter  and  darker  rings  which  are  larger  in  the  center  and 

become  narrower  as  the  off-axis  distance  increases.  This  occurs  in  every  frame,  no 

matter  what  the  scene;  however  it  becomes  most  apparent  when  viewing  large  highly 
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uniform  regions.  It  should  be  noted  that  the  k-means  algorithm  defines  end  members 
based  on  both  spectral  correlation  and  intensity  differences.  So  while  there  are  only  four 
true  spectral  end  members  in  the  scene,  regions  such  as  the  sky  can  be  further  broken  into 
groups  by  intensity.  This  feature  is  also  what  allows  the  algorithm  to  see  the  fringe 
pattern.  Spectrally,  the  two  fringing  end  members  are  the  same,  but  their  intensity  is 
slightly  different.  In  fact,  it  is  hard  to  distinguish  the  two  fringe  intensities  on  the  bottom 
plot  because  they  are  almost  right  on  top  of  each  other  with  less  than  a  0.5  K  difference. 
This  is  close  to  the  noise  equivalent  delta  temperature  (NEDT)  for  the  instrument  and 
therefore  should  not  be  visible  in  most  cases.  Other  end-member  algorithms  exist  that 
perform  grouping  based  solely  on  spectral  features.  Using  this  type  of  algorithm  would 
eliminate  the  subdivision  of  these  spectrally  uniform  regions. 
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Figure  17.  Results  of  the  end  member  analysis  after  the  spectral  axis  scaling  has  been  corrected. 
Now,  all  spectral  features  line  up  correctly  and  there  are  no  concentric  rings  due  to  off-axis  scaling. 
What  is  visible,  is  the  fringe  pattern  created  by  the  interferometer  itself.  This  is  an  artifact  of  forcing 
10  end  members  when  there  are  only  four  true  end-members.  The  light  and  dark  rings  are  separated 
by  a  temperature  difference  of  less  than  0.5  K  which  is  close  to  the  NEDT  of  the  instrument. 

Ten  end  members  were  used  above  to  accentuate  the  off-axis  scaling  issues. 

Figure  18  and  Figure  19  show  the  same  analysis  done  with  a  more  realistic  five  end 
members.  Figure  18  shows  the  end-member  analysis  before  spectral  axis  correction. 
Again,  the  concentric  rings  in  the  sky  region  are  visible  and  in  this  case  the  road  is 
incorrectly  grouped  with  another  end-member.  Figure  19  shows  the  end-member 
analysis  after  the  spectral  axis  correction.  Now  the  sky  is  correctly  grouped  by  altitude 
and  the  road  is  identified  as  its  own  end-member.  This  illustrates  the  importance  of  the 
spectral  axis  correction  for  applications  such  as  anomaly  detection  or  target  detection. 

For  example,  if  an  intelligence  analyst  were  using  algorithms  such  as  these  to  search  for 
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an  anomaly  in  the  scene  and  they  didn't  apply  a  proper  off-axis  spectral  calibration,  then  a 


key  feature  in  the  scene  might  be  missed. 


Figure  18.  End  member  analysis  before  spectral  axis  correction  using  five  end  members.  The 
concentric  rings  are  still  visible  in  the  sky  region  and  the  road  is  incorrectly  grouped  together  with 

another  end-member. 


Figure  19.  End  member  analysis  after  spectral  axis  correction  using  five  end-members.  The  sky  is 
now  appropriately  grouped  by  altitude  and  the  road  is  properly  identified  as  its  own  end  member. 
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IV.  Polarimetric-HSI  Calibration 


In  addition  to  radiometric  and  off-axis  spectral  calibration,  a  polarimetric 
calibration  is  necessary  in  order  to  fully  exploit  the  spectro-polarimetric  dataset.  This 
necessity  arises  from  the  fact  that  the  gain  of  the  system  is  sensitive  to  the  instrument 
polarizer  angle.  Two  factors  must  be  considered:  first,  the  instrument  polarizer  is  self- 
emitting  polarized  light  in  the  LWIR.  This  distorts  the  polarization  state  of  light  from  the 
scene  that  strikes  the  detector  since  the  emitted  light  is  orthogonally  polarized  with 
respect  to  the  transmitted  polarized  light  from  the  scene.  Second,  it  is  possible  that  the 
other  optics  in  the  system  have  a  slight  polarization  preference,  although  this  should  be 
small  for  transmissive  optics.  The  beam  splitter  in  the  interferometer  is  likely  to  have  the 
greatest  polarization  effect  due  to  its  higher  reflectance  and  angle  of  incidence  with 
respect  to  incoming  light.  To  achieve  highly  accurate  Stokes  vectors,  the  full  spectro- 
polarimetric  gain  of  the  system  must  be  determined  for  each  pixel. 

The  first  part  of  this  chapter  starts  by  developing  a  mathematical  framework  for 
use  as  a  polarimetric  system  model  using  Stokes  vectors  and  the  Mueller  calculus.  This 
polarimetric  system  model  is  then  used  to  derive  the  calibration  methodology  based  on 
the  current  assumption  that  the  instrument  polarizer  is  ideal,  or  in  other  words  its 
diattenuation  as  described  by  Equation  (21)  is  equal  to  one.  The  current  polarimetric 
calibration  procedure  involves  performing  a  radiometric  calibration  at  each  instrument 
polarizer  angle.  It  will  be  shown  that  under  the  ideal  polarizer  assumption,  this 
polarimetric  calibration  procedure  works  perfectly  well. 


45 


The  next  section  will  then  discuss  the  polarimetric  calibration  without  the 
assumption  of  an  ideal  polarizer.  In  this  section,  the  experimental  methods  required  for 
characterizing  the  non-ideal  nature  of  the  system  and  providing  a  calibration  for  these 
non-idealities  will  be  demonstrated. 

Next,  the  Analysis  of  results  section  quantifies  how  ideally  the  system  is 
performing  and  compares  the  accuracy  of  the  measurements  using  the  current  calibration 
method  versus  the  accuracy  of  the  measurements  with  the  additional  correction  for  a  non¬ 
ideal  polarizer.  It  also  includes  a  discussion  of  the  noise  and  error  in  the  measurements. 

Next,  the  NESR  of  the  system  is  determined  using  an  unpolarized  blackbody 
source.  The  concept  of  noise  equivalent  spectral  DoLP  (NESDoLP)  is  then  discussed  and 
experimentally  determined  using  the  same  data  set.  These  same  parameters  are  then 
determined  for  a  polarized  source  and  the  results  from  both  calibration  methods  are 
compared. 

Polarimetric  System  Model 

To  describe  the  system,  a  Mueller/Stokes  representation  is  used  where  each 
optical  component  in  the  system  can  be  described  by  its  own  Mueller  matrix  and  offset 
Stokes  vector.  The  first  optical  component  that  the  scene  radiance  will  encounter  is  the 
polarizer.  The  Mueller  matrix  of  the  polarizer  will  operate  on  the  scene  Stokes  vector, 
creating  a  new  Stokes  vector.  This  new  Stokes  vector  is  added  to  the  polarizer's  offset 
Stokes  vector.  It  will  then  encounter  the  next  optical  element  in  the  system  which  will 
then  convert  that  Stokes  vector  into  another  Stokes  vector  and  so  on  until  the  final  Stokes 
vector  reaches  the  detector.  A  simplified  schematic  of  the  system  is  shown  in  Figure  20 
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where  an  input  Stokes  vector  from  the  scene  travels  through  the  system  getting  modified 
by  each  optical  component  in  the  system.  After  passing  through  the  system,  a  Stokes 
vector  reaches  the  detector  which  is  a  function  of  the  system  Mueller  matrix  and  the 
instrument  offset  Stokes  vector. 


Scene  Stokes  Vector 


M,  (6p  ) 


Detector 

S  d  (®p) 


Figure  20.  Simplified  P-HSI  system  where  an  input  Stokes  vector  from  the  scene  travels  through  the 
system  getting  modified  by  each  optical  component  in  the  system.  After  passing  through  the  system, 
a  Stokes  vector  reaches  the  detector  which  is  a  function  of  the  system  Mueller  matrix  and  the 

instrument  offset  Stokes  vector. 


Using  this  approach,  the  Stokes  vector  reaching  the  detector  after  having  passed 
through  each  optical  component  is  given  by 

S/»,)=M„(...(M2(M,(e,)S,+S1(ep))+S2)+...)  +  S.  (38) 

where  Mj  (0  )  is  the  Mueller  matrix  of  the  polarizer,  Ss  is  the  Stokes  vector  of  the 
incoming  scene  radiance,  Mn  is  Mueller  matrix  of  the  n,h  optical  component,  and  Sn  is 
the  Stokes  vector  containing  the  self-emission  of  the  nth  optical  component.  Ss  can  be 
written  as 


S.  = 


so  ~  1 


V  S3 


■L 


(39) 
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where  s0,  .v, ,  s2,  and  s3  are  the  normalized  Stokes  parameters  and  Ls  is  the  scene 
radiance.  Notice  that  Mj(#  )is  the  Mueller  matrix  of  the  polarizer  and  is  a  function  of 

the  polarizer  angle,  6  p ,  and  5,  (6p)  is  the  self-emission  of  the  polarizer  which  is  also  a 

function  of  the  polarizer  angle.  It  should  also  be  noted  that  each  term  in  Equation  (38)  is 
a  spectral  quantity  and  therefore  is  wavenumber  dependent;  however,  the  notation  is 
dropped  for  brevity.  Equation  (38)  captures  the  cumulative  transmitted  and  self-emitted 
radiances  of  each  optical  component  as  well  as  the  polarization  state  as  light  travels 
through  the  optical  system.  If  Equation  (38)  is  multiplied  out,  it  can  be  written  as 

S/0p)  =  M(0p)S,+S„(fl,)  (40) 

where  M(6U)  is  the  Mueller  matrix  of  the  system  as  a  whole  and  So(6p)  is  the  cumulative 
system  offset  Stokes  vector.  Utilizing  product-sum  notation,  M (0  )  can  be  written  as 

(41) 

v — ! - v - j  Polarizer 

HSI 

where  the  first  term  is  the  Mueller  matrix  of  the  HSI  system  (not  including  the  polarizer) 
and  the  second  term  is  the  Mueller  matrix  of  the  polarizer.  Equations  (38)  and  (40) 
describe  the  Stokes  vector  arriving  at  the  detector.  Upon  striking  the  detector  the  Stokes 
vector  is  converted  from  radiance  units  into  arbitrary  digital  units  (DN).  The  detector  is 
also  insensitive  to  the  polarization  state  of  incident  light  and  therefore  only  measures  the 

s0  element  of  Sd  (0  ) .  To  account  for  this,  a  detector  vector  can  be  defined  as 

D  =  (d  0  0  0) 

where  d  is  proportional  to  the  spectral  responsivity  of  the  detector. 
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(42) 


The  measured  signal  is  then  given  by 

YmWr)  =  D-Sd(6p)  =  b-M(ep)S,+Y^0p)  [DN-cm-1]  (43) 

where  Ya  (9  )  is  the  offset  in  DN  units.  Plugging  in  for  D,  M(#p),  and  Ss  gives 

™oo  (Op)  moi  (Op)  m<pSe„)  nh:S°py 
"ho (Op)  mil  (Op)  ml2(0p)  ™l3(0p) 

/Wjo (6 p )  fn2i  {Op)  n?22  ( 0 p )  ^23 (^>) 

m30(^p)  m3|C^)  '«32C^)  m33(^)y 

which  gives, 

y«  (#„ )  =  d  (™oo  (#P )  +  ^1 W01  (#p )  +  s2m01 (1 9p )  +  s3m03  (0p  j)Ls+Y0  {0p  ) 
and  factoring  out  m00  gives, 

r,„  (9t)  =  d-  )(l  +  4,1%  (g, )  +  a,  rffe  (g„ )  +  (0,  ))L,  +  y,  (g,)  (44) 

Gain  Offset 

where  m01 ,  m02 ,  m03  are  the  normalized  Mueller  elements  of  the  system  as  a  whole. 

Equation  (44)  describes  the  measured  signal  at  the  detector  as  a  function  of 
polarizer  angle  and  the  first  row  of  the  Mueller  matrix  of  the  system.  As  described  in 
Chapter  III,  radiometric  calibration  involves  inputting  unpolarized  blackbody  radiation 
into  the  system  at  two  different  temperatures  and  solving  for  the  gain  and  offset.  When 
this  is  done,  s{,  s2,  and  s3  are  zero  and  Equation  (44)  reduces  to  Equation  (23)  where  the 

gain  is  now  written  in  terms  of  the  spectral  detector  responsivity,  described  with  d ,  and 
the  spectral  transmission  of  the  optical  system,  described  with  mm  (6 p ) .  After  solving  for 

the  gain  and  offset,  as  discussed  in  Chapter  III,  and  converting  to  calibrated  radiance, 
Equation  (44)  can  be  written  as 

Lm  (0P  )  =  (l  +  s1th01  (0p  )  +  s2mm  (9p  )  +  s3m03  (9p  ))Ls 
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Ym(0p)  =  (d  0  0  0). 


(45) 


where  Lm  (6p )  is  the  calibrated  radiance  measured  by  the  detector.  Notice  that  there  are 

four  unknown  Stokes  parameters  from  the  scene  that  must  be  determined.  This  requires 
measurement  of  the  scene  at  four  different  polarizer  angles.  When  this  is  done,  the 


resulting  system  of  equations  based  on  Equation  (45)  can  be  written  in  matrix  form  as 

'  4(0)  ^  fl  m01( 0)  m02 (0)  ^(0)  UP 

4,(45)  =  1  m01(45)  m02  (45)  m03(45)  y  ^ 

4(90)  1  m01  (90)  m02  (90)  m03(90)  '  s2  * 

V4(135)J  [l  m01(135)  m02(135)  m03(135)J 


or 

4=W  -S,.  (47) 

In  Equation  (47),  Lm  is  called  the  channel  vector  and  it  contains  the  measured  radiance 
at  each  of  the  four  polarizer  angles  (or  channels),  W  is  called  the  system  matrix,  and  Ss 

is  the  input  Stokes  vector  from  the  scene.  Clearly,  the  system  matrix  transforms  the  input 
Stokes  vector  into  the  measured  radiance  at  each  of  the  four  polarimetric  channels.  The 
system  matrix,  therefore,  represents  the  polarimetric  effect  that  the  system  has  on  the 
input  signal.  The  goal  of  the  polarimetric  calibration  is  to  determine  the  system  matrix  so 
that  its  effects  can  be  removed.  Once  the  system  matrix  is  known,  the  input  Stokes  vector 
can  be  determined  by  multiplying  both  sides  of  the  equation  by  the  pseudo-inverse  of  the 
system  matrix,  that  is 

S,  =  w+4  (48) 

where  W+  is  the  pseudo-inverse  of  W .  In  the  literature,  W+  is  known  as  the  data 
reduction  matrix,  R  so  that 
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Ss=*-Lm. 


(49) 


Up  until  this  point  no  assumptions  have  been  made  about  the  polarizer.  The  next 
section  will  determine  the  system  matrix  under  the  assumption  that  the  instrument 
polarizer  is  ideal.  After  that,  the  more  general  case  of  a  non-ideal  polarizer  is  developed. 


Polarimetric  Calibration  under  the  Ideal  Polarizer  Assumption 

To  determine  the  system  matrix,  the  form  of  mm(0p) ,  m()[ (&p),  m02 (&p),  and 

mO3(0p) must  be  determined.  This  can  be  done  by  carrying  out  the  multiplication  in 


Equation  (41).  For  an  ideal  polarizer,  Mj(0  )  is  given  by 


f  1 


**>(*„)  = 


1 

2 


cos(26»p) 

sin(26»p) 


cos(2  0p) 
cos2(26»p) 
sin(2#p)cos(2#p) 
0 


sm(20p)  0^ 

sin(2#p)cos(2#p)  0 
sin2  (2  0p)  0 

0  0y 
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Plugging  this  into  Equation  (41)  gives 


M(^)  =  M^M1(^)  = 


1 

K 

m02 

m03 

1 

cos(26»p) 

sin(2  9p) 

(T 

K o 

m[  i 

mh 

m'l3 

cos(2  ep) 

cos2  (20 p ) 

sin(2d;Jcos(2d;J 

0 

1 

m20 

m2\ 

m22 

m23 

sin  (20p) 

sin(2d;) )  cos(2d;j ) 

sin2(2/9p) 

0 

2 

m30 

K 

m32 

m33, 

o 

0 

0 

0, 

where  M^sv  is  the  normalized  Mueller  matrix  of  the  HSI  system  not  including  the 
polarizer  and  M,  (6p )  is  the  Mueller  matrix  of  the  polarizer.  Carrying  out  this 
multiplication  for  the  first  row  only  and  normalizing  by  the  m00  element  gives 


M (0P )  =  \ (l  +  K\  cos(20p )  +  m'02  sin(2 0p)) 


1 

cos(2  0p) 

sin(26»p) 

0^ 

>KSeP) 

0 

™20  (Op) 

™oJQp) 

m22(dp) 

0 

"bo  (Op) 

K(0P) 

KiVp) 
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Therefore,  m00 (0p)  =  ^(l  +  m'm  cos(26>p )  +  m'02  s\n(20p )),  m(n (0p )  =  cos(20p )  , 
m02  {6 p )  =  sin(2#p ) ,  and  m03  ($/;)  =  0.  Notice,  that  since  m(X,  (0p )  =  0 ,  ,s’3  cannot  be 


measured  with  this  system.  This  result  makes  sense  for  a  linear  polarizer  because  s3 


represents  circularly  polarized  light.  Plugging  these  into  the  system  matrix  gives 


W 

ideal 


11  0^ 
1  0  1 

1  -1  0 


1  0  -1 


(51) 


Solving  for  the  data  reduction  matrix  by  taking  the  pseudo-inverse  of  Wideal  gives 


R 


ideal 


^0.25  0.25  0.25  0.25^ 
0.5  0  -0.5  0 

v  0  0.5  0  -0.5, 


The  input  Stokes  vector  can  then  be  determined  according  to  Equation  (49)  as 


(52) 


.v. 


^0.25  0.25  0.25  0.25 
0.5  0  -0.5  0 

,0  0.5  0  -0.5 


A 

f  4(0)  ^ 

4(45) 

4(90) 

) 

14  (135)  J 

(53) 


Carrying  out  this  multiplication  gives 


Lm  (0)  +  Lm  (45)  +  Lm  (90)  +  Lm  (135) 


(54) 


4(0) -4  (90 ) 


(55) 


and 


_  4  (45)  ~  4  (1 35) 


(56) 
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This  result  indicates  that  if  the  polarizer  is  ideal,  then  polarimetric  calibration  can 
be  achieved  simply  by  performing  a  radiometric  calibration  at  each  of  the  four  polarizer 
angles.  These  calibrated  radiances  can  then  be  plugged  into  the  above  equations  to 
determine  the  Stokes  parameters  of  the  scene. 

If  the  polarizer  is  not  perfectly  ideal,  however,  the  elements  of  the  system  matrix 
will  not  be  given  by  simple  trigonometric  functions  and  some  residual  error  in  the 
polarimetric  measurements  will  remain  if  the  above  equations  are  used.  The  magnitude 
of  the  residual  error  will  depend  on  how  non-ideal  the  polarizer  is  and  how  large  of  an 
effect  the  optical  components  in  the  HSI  system  have  on  the  polarization  of  light 
travelling  through  it.  For  a  non-ideal  polarizer,  the  system  matrix  must  be  experimentally 
determined.  The  next  few  sections  walk  through  the  experimental  determination  of  the 
system  matrix  and  the  analysis  of  the  system  performance. 

Polarimetric  Calibration  Without  the  Ideal  Polarizer  Assumption 

In  order  to  solve  for  the  system  matrix  in  Equation  (47),  the  channel  vector  must 
be  measured  for  multiple  input  Stokes  vectors.  The  number  of  input  Stokes  vectors  must 
be  large  enough  to  allow  for  a  proper  least  squares  solution  to  the  resulting  matrix 
equation.  In  this  case,  Equation  (47)  becomes  the  following  matrix  equation: 

Lm=W-S  (57) 

where  the  channel  vector  is  now  a  M  x  N  matrix,  Lm ,  where  M  is  the  number  of 

instrument  polarizer  angles  and  N  is  the  number  of  unique  Stokes  input  vectors.  It 

follows  that  S  is  now  a  4xN  matrix  where  each  row  corresponds  to  an  element  in  the 

Stokes  vector.  In  this  work,  M=4,  which  corresponds  to  measurements  at  0,  45°,  90°,  and 
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135°.  Clearly,  in  order  to  achieve  an  accurate  solution  for  W ,  the  input  Stokes  vectors 
must  be  known  with  a  high  accuracy.  This  presents  a  challenge  when  working  in  the 
LWIR  for  reasons  which  will  be  discussed  below.  An  experimental  method  to  overcome 
this  challenge  and  generate  highly  accurate  input  Stokes  vectors  was  originally  developed 
by  Persons,  C.  M.  et  a.l  in  2007  (see  Reference  [20]).  The  methodology  used  below  is 
based  on  the  method  used  in  that  paper. 

Methodology  Overview 

An  external  polarizer  is  used  to  generate  the  known  input  Stokes  vector;  however, 
in  the  LWIR,  knowledge  of  the  input  Stokes  vector  is  more  complicated  than  in  the 
visible  region  due  to  the  self-emission  and  reflected  components  of  the  external  polarizer. 
Figure  21  illustrates  this  concept  where  S£is  the  self-emitted  component,  S7  is  the 

transmitted  component,  and  SR  is  the  reflected  component. 


Figure  21.  Experimental  setup  showing  blackbody  radiation  passing  through  a  linear  polarizer  and 
reaching  the  Telops  Hypercam.  Illustrates  the  self-emitted  and  reflected  components  that  also  reach 
the  instrument  making  it  difficult  to  determine  the  transmitted  component. 

In  order  to  overcome  this  difficulty,  the  self-emitted  and  reflected  components 
must  be  effectively  removed  to  isolate  the  transmitted  component.  If  they  are  not 
removed,  they  will  negatively  affect  the  data  in  two  ways.  First,  the  amount  of  radiance 
being  emitted  and  reflected  from  the  polarizer  is  likely  not  known  with  great  accuracy 
and  therefore  knowledge  of  the  input  Stokes  vector  would  be  inaccurate.  This  would 
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introduce  error  in  the  calibration.  Second,  the  reflected  and  emitted  radiance  from  the 


polarizer  are  orthogonally  polarized  with  respect  to  the  transmitted  component.  This  has 
the  effect  of  reducing  the  DoLP  of  the  signal  that  makes  it  to  the  instrument  and  therefore 
reducing  the  SNR  of  the  polarimetric  measurements.  To  remove  the  reflected  and 
emitted  components,  the  setup  shown  is  modified  with  an  additional  reference  source  to 
control  the  reflected  component  and  data  is  then  collected  at  two  different  target  source 
temperatures  so  that  the  two  can  be  differenced.  This  reference  source  overfills  the 
polarizer  to  ensure  a  uniform  reflected  component.  Figure  22  illustrates  the  modified 
setup.  In  this  setup,  the  external  polarizer's  plane  is  rotated  by  4°  and  the  reference 
blackbody  is  placed  at  the  specular  angle  of  8°  with  respect  to  the  Telops  instrument.  The 
reference  blackbody  is  then  kept  at  a  constant  temperature  in  order  to  control  the  reflected 
component.  During  the  experiment,  it  was  found  that  a  baffle  that  prevented  light  from 
the  target  source  from  reflecting  off  of  the  metallic  edge  of  the  reference  blackbody  and 
then  off  of  the  external  polarizer  into  the  system  was  also  necessary. 


Figure  22.  Experimental  setup  used  for  polarimetric  calibration.  Two  different  target  source 
temperatures  are  used  for  differencing  and  the  reference  source  is  placed  at  the  specular  angle  and 

kept  at  a  constant  temperature. 
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With  this  setup  in  place,  the  external  polarizer  is  used  to  generate  the  input 
polarization  states,  and  at  each  polarization  state,  data  is  collected  at  two  different 
temperatures.  The  input  Stokes  vectors  are  then  differenced  so  that 

S2  =  ST !  +  pP  +  pP 

-  S,  =STi+^  +  X  .  (58) 

AS  =  A  ST 

This  gets  rid  of  the  emitted  and  reflected  components  and  isolates  the  transmitted 
component.  After  differencing,  the  Stokes  vector  that  reaches  the  instrument  is 
simplified  to 

1 

Z)cos(26f) 

D  sin(2$) 

0 

where  D  is  the  diattenuation  of  the  external  polarizer  and  A Ls  is  the  differenced 
radiance  of  the  blackbody  after  transmission  through  the  external  polarizer.  Note  that 
since  s3  =  0 ,  it  will  be  dropped  and  the  three-element  Stokes  vector  will  be  used  from 
now  on.  Now,  Equation  (57)  becomes 

ALm  =  W  •  AS  .  (60) 

For  example,  if  only  four  input  polarization  states  were  used  where  the  external  polarizer 
was  oriented  at  0°,  45°,  90°,  and  135°,  then  the  resulting  equation  would  be 
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(61) 


where  ALm  (#p )  is  the  measured  radiance  after  calibration  and  after  differencing,  and  the 

superscript  corresponds  to  the  external  polarizer  angle  (and  therefore  the  input  Stokes 
vector).  Notice  that  A Ls  must  be  determined  before  W  can  be  determined.  To 

determine  A Ls  the  channel  vectors  in  ALm  are  averaged  and  the  Stokes  vectors  in  AS 
are  averaged  to  generate  the  following  equation: 


-( AL  0  +  AL  45  +  AL  90  +  AL  135 )  =  W • 

^  y  m  m  m  m  J 


o 

v  0  y 


(62) 


In  the  derivation  of  the  system  matrix,  it  was  shown  that  the  first  column  of  W  is  a 
column  of  ones.  This  was  simply  a  result  of  performing  a  radiometric  calibration  at  each 
of  the  four  instrument  polarizer  angles.  Utilizing  this,  the  right  side  of  Equation  (62)  can 
be  multiplied  out  to  give 


-(AL  0  +  AL  45  +  AL  90  +AL  135  )  = 

a  y  m  m  m  m  J 


^AL  ^ 

AL. 
A Ls 


(63) 


which  means  that  the  best  statistical  solution  for  A Ls  is  just  the  average  of  the  channel 
matrix,  or 

AL  =  T£(AZ?+AZ;4J+AZr  +  AZ;35).  («) 
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Now,  everything  that  is  necessary  to  compute  the  system  matrix  is  available.  Solving  for 


W  gives 


W  =  AL  -AS4 


(65) 


or 


W 


'  <0 0) 

<(0) 

<(0) 

<5(0)  ' 

<(45) 

<(45) 

<(45) 

<5(45) 

<(90) 

<(90) 

<(90) 

ALf  (90) 

.<(135) 

<(135) 

<(135) 

<T  (135), 

A U 


Oil  1  ^ 
D  0  -D  0 
0  D  0  -D 


(66) 


where  again,  AS4  is  the  pseudo-inverse  of  AS  .  From  this,  the  data  reduction  matrix,  R  , 
can  be  determined  as  before  and  this  can  be  used  to  convert  the  measured  channel  vector 


of  any  scene  to  the  corresponding  input  Stokes  vector  without  the  assumption  of  an  ideal 
instrument  polarizer. 

In  the  example  above,  only  four  input  Stokes  vectors  were  used  to  make  the  math 
conceptually  simpler.  By  using  additional  input  Stokes  vectors  that  are  mixtures  of  the 
four  primary  ones  described  above,  the  calibration  should  be  more  accurate  for  a  greater 
number  of  input  states.  In  light  of  this,  12  input  Stokes  vectors  were  used  rather  than  4 
by  incrementing  the  external  polarizer  by  15°  starting  at  0°  and  going  to  165°. 

Experimental  Setup 

The  physical  setup  used  for  this  experiment  was  briefly  discussed  above  and  is 
shown  in  Figure  22.  This  section  describes  the  experimental  setup  in  much  greater  detail. 
The  two  target  source  temperatures  were  55  °C  and  65  °C  and  the  reference  blackbody 
was  kept  at  30  °C.  The  2-inch  target  blackbodies  have  a  temperature  stability  of  0.003  °C 
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and  the  stability  of  the  6-inch  reference  blackbody  is  0.5  °C.  The  onboard  blackbodies 
used  for  calibration  were  kept  at  20  °C  and  60  °C  and  the  integration  time  was  kept  at 
200//.s  for  all  collections.  At  these  temperatures,  the  calibration  blackbodies  bounded  the 
scene  radiance  and  stayed  within  the  linear  range  of  the  instrument.  The  data  was 
collected  at  16  cm-1  spectral  resolution  and  64  acquisitions  of  each  scene  were  collected, 
as  well  as  34  acquisitions  for  each  calibration  blackbody  collection  for  averaging 
purposes.  The  Hypercam  FPA  was  windowed  down  so  that  only  the  polarizer  and  part  of 
the  target  blackbody  was  in  the  scene.  The  broadband  image  of  the  scene  is  shown  in 
Figure  23.  As  mentioned  before,  a  baffle  was  placed  on  the  left  side  of  the  polarizer  so 
that  light  from  the  target  blackbody  would  not  reflect  off  of  the  reference  blackbody  and 
then  off  of  the  polarizer  into  the  instrument.  Before  the  baffle  was  put  into  place,  a  stripe 
could  be  seen  in  the  image  of  the  polarizer  where  the  reflection  off  of  the  polarizer  was 
taking  place.  The  experimental  setup  parameters  are  outlined  in  Table  2.  Parameters  that 
are  colored  red  are  thought  to  be  the  primary  limitations  of  the  experimental  setup  which 
could  be  improved  upon  in  future  experiments. 


Table  2.  Experimental  Setup  Parameters 


Parameter 


Specification 


Target  source  temperatures 
Target  source  stability 
Reference  source  temperature 
Reference  source  stability 

Calibration  source  temperatures 

Avg.  external  polarizer 
diattenuation 


55,  65  °C 
0.003  °C 
30  °C 
0.5  °C 
20,  60  °C 


0.992 
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Integration  time 

Spectral  range 

Spectral  resolution 

Polarizer  type 

Polarizer  material 

Polarizer  orientation  accuracy 

Frame  averaging 


200  jus 
900-1250  cm'1 
16  cm'1 


Wire-grid  (linear) 
ZnSe 


1° 
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In  order  to  ensure  proper  alignment  of  the  external  and  instrument  polarizers,  a 
dataset  was  collected  where  the  polarizer  was  held  constant  and  the  external  polarizer 
was  rotated  from  0  to  180°  in  4°  increments  except  near  the  maximum  intensity  point, 
where  it  was  rotated  in  1°  increments  for  higher  resolution.  This  was  done  to  map  the 
angle  of  the  rotation  stage  holding  the  polarizer  to  the  transmission  axis  of  the  polarizer 
itself.  The  results  are  shown  in  Figure  24.  Since  the  instrument  polarizer  was  kept  in  the 
horizontal  position,  the  angle  of  the  external  polarizer  at  which  maximum  signal  was 
achieved  was  used  to  determine  the  proper  horizontal  alignment  of  the  external  polarizer. 
A  cos2  (0)  curve  was  fit  to  the  data  using  Malus'  law  which  is  given  by 


4,  =4  cos2  (0). 


(67) 


This  allowed  for  proper  orientation  of  the  external  polarizer  for  generating  the  necessary 
input  Stokes  vectors. 
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Rotation 


Polarizer  with  blackbody 
behind  it 


Figure  23.  Band-averaged  image  of  the  polarizer  with  the  65  °C  blackbody  in  the  background.  Note 
the  active  area  of  the  blackbody  overfills  the  polarizer  but  not  the  entire  rotation  stage.  The  baffle 

can  be  seen  on  the  left  side  of  the  image. 


0  20  40  60  80  100  120  140  160  180 

External  Polarizer  Orientation  [Deg] 


Figure  24.  Data  collected  from  the  external  polarizer  as  it  was  rotated  through  180°.  The  angle  of 
the  external  polarizer  at  which  maximum  signal  was  received  corresponds  to  the  best  alignment 
between  the  external  and  internal  polarizer.  The  solid  line  is  the  best  fit  cos2  curve. 


Using  this  setup,  data  was  collected  in  15°  increments  starting  at  0°  and  going  to 

165°  degrees  resulting  in  12  input  Stokes  vectors  at  the  cold  temperature  and  12  input 

Stokes  vectors  at  the  hot  temperature.  It  was  important  to  take  the  high-temperature 
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measurement  immediately  following  the  low  temperature  measurement  in  order  to 
minimize  the  temperature  drift  of  the  external  polarizer.  Therefore,  two  target 
blackbodies  were  used  with  one  set  at  the  high  temperature  and  one  set  at  the  low 
temperature  so  they  could  be  switched  out  quickly.  The  sequence  of  instrument  polarizer 
angles  used  was  also  optimized  in  order  to  minimize  temperature  drift  between 
orthogonal  polarizer  measurements,  i.e.,  the  instrument  polarizer  collection  order  was  0°, 
90°,  45°,  and  then  135°.  This  method  is  implemented  due  to  the  fact  that  calculation  of 
the  .s',  and  .v2  Stokes  elements  requires  differencing  the  two  orthogonal  states.  Figure  25 

shows  the  s0,  .s, ,  and  ,s2  Stokes  images  from  the  scene  with  the  external  polarizer  at  15° 
and  the  target  blackbody  at  65  °C.  In  the  ,v0  image,  the  background  scene  including  the 
blackbody  and  the  baffle  is  present,  whereas  in  the  .s,  and  .v2  images  the  background 
scene  is  black.  This  indicates  that  there  was  very  little  temperature  drift  in  the  scene 
during  the  collection.  Additionally,  Figure  26  shows  the  ,s0  image  after  differencing  the 
hot  and  cold  temperatures.  This  shows  that  the  background  scene  is  no  longer  visible 
even  in  the  ,v0  image,  which  indicates  that  there  was  minimal  temperature  drift  between 


collection  of  the  hot  and  cold  datasets. 


Figure  25.  Stokes  images  of  the  scene  with  the  target  blackbody  at  65  °C.  The  background  scene  is 
visible  in  the  s«  Stokes  image,  whereas  it  is  not  visible  in  the  Si  and  s2  Stokes  images  indicating  that 
there  was  very  little  thermal  drift  in  the  scene  during  the  collection. 


62 


Figure  26.  s0  image  after  differencing  the  hot  and  cold  data.  The  dark  background  scene  indicates 
that  the  thermal  drift  between  collection  of  the  hot  and  cold  data  was  small. 

Next,  a  channel  matrix  was  generated  based  on  the  calibrated  differenced  data. 
Figure  27  provides  a  useful  way  of  visualizing  the  channel  matrix  which  helps  to 
understand  its  physical  significance.  In  the  figure,  each  column  is  a  channel  vector  which 
corresponds  to  an  input  Stokes  vector.  For  visualization  purposes,  only  four  input  Stokes 
orientations  are  shown  rather  than  all  12. 
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Figure  27.  Visualization  of  the  channel  matrix  where  each  column  is  a  channel  vector  corresponding 

to  an  input  Stokes  vector. 


In  order  to  use  the  data  for  calculation  of  the  system  and  data  reduction  matrices, 
only  a  circular  region  within  the  polarizer  was  used.  To  do  this,  a  circular  mask  was 
created  to  remove  the  scene  other  than  the  circular  region  within  the  external  polarizer. 
This  is  shown  in  Figure  28. 


Figure  28.  Circular  mask  to  remove  all  scene  data  other  than  the  circular  region  within  the 

polarizer. 
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With  this  mask  in  place,  it  is  easy  to  show  the  effect  that  differencing  the  hot  and  cold 
data  has  on  .s', ,  s2,  and  DoLP.  Figure  29  shows  the  .v, ,  s2,  and  DoLP  images  of  the 
polarizer  when  it  is  horizontally  aligned  with  the  instrument  polarizer.  This  should  result 
in  an  sr  close  to  one,  an  s2  close  to  zero,  and  a  DoLP  close  to  100%.  This  proves  to  be 
true  for  the  differenced  data,  while  it  is  far  from  true  for  the  undifferenced  data  due  to  the 
reflected  and  emitted  component  discussed  before.  For  the  undifferenced  data,  .v,  is 
approximately  0.2  and  the  percent  DoLP  is  approximately  20%.  For  the  differenced  data, 
s1  is  approximately  0.995  and  the  percent  DoLP  is  approximately  99.5%.  In  practice,  a 

DoLP  of  100%  is  not  achievable  since  there  will  always  be  some  diattenuation  that  drives 
it  down;  however,  noise  in  the  measurements  can  potentially  cause  a  DoLP  of  over  100% 
when  using  the  Modified  Pickering  method.  To  demonstrate  this  mathematically, 
consider  the  following  calculation  of  s1  that  includes  noise: 


s,  = 


(Lm  (0)  +  n(0))  ~(Lm  (90)  +  «(90)) 


\[(Lm  (0)  +  n(  0))  +  (Lm  (90)  +  n(90))  +  (Lm  (45)  +  n(45))  +  (Lm  (135)  +  n(135))] 


(Lm  (0)  -  Lm  (90))  +  (n(0)  -  n(90)) 


.  (68) 


\[{Lm  (0)  +  Lm  (90)  +  Lm  (45)  +  Lm  (135))  +  (n(  0)  +  n(90)  +  n(45)  +  n(135))] 


Therefore,  if  Lm  (90)  is  close  to  zero  as  it  should  be  if  the  external  polarizer  is  horizontal, 
and  if  [ Lm  (45)  +  Lm  (135)]  <  [ Lm  (0)  +  Lm  (90)]  due  to  noise,  then  this  could  result  in  an 


artificially  high  calculation  of  .v, .  This  noise  is  minimized  by  averaging  multiple  frames, 

but  in  this  case,  thermal  drift  of  the  external  polarizer  between  measurement  of  the  hot 
and  cold  data  sets  can  have  a  similar  effect.  It  is  also  interesting  and  somewhat 
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counterintuitive  to  note  that  although  the  s1  in  the  undifferenced  data  is  much  smaller 
than  in  the  differenced  data,  the  s2  value  in  the  undifferenced  data  is  actually  closer  to 
zero  than  in  the  differenced  data.  Another  look  at  the  math  will  reveal  that  this  is  a 
perfectly  valid  result.  For  the  undifferenced  data,  s1  is  given  by 

_  (>? (0)r  +  5(0) ^  +  5(0)£ ) - (5 (90)r  +  5 (90) ^  +  5 (90) E ) 

'V|  ~  (s(o)T+s(o)R+s(o)E)+(smT+smR+smE) ' 

Assuming  5(90)y  «  0  ,  5(0)  «  0,  and  S(0)E  «  0  this  gives 


s,  _'S(Q)T-smR-smE 
Sl~  s«f)T+smR+smE 

whereas  with  the  difference  data,  s1  is  given  by 


(69) 


AS(0)r 

AS(0)r 


(70) 


This  shows  that  in  the  differenced  data,  s1  should  be  close  to  1,  but  in  the  undifferenced 


data  if  the  emitted  and  reflected  components  are  significant  then  s1  will  not  be  close  to 


one.  Calculation  of  s2  in  the  undifferenced  data  is  given  by 

(5(45)r  +S(45)r  +  S(45)£)-(5;(135)r  +5(135)*  +5(135)£) 
2  (S(45)r  +  S(45)r  +  5(  45)  £ )  +  (5(135)T  +  5(135)  R  +  5(135)  E ) 


where  5(45)T  *  5,(135)T ,  5,(45)ff  «  5(135)ff  ,  and  5( 45) £  «  5(135) £ ,  which  gives 


s 


2 


0 

2(St  +Sr+Se ) 


«  0 . 


For  the  differenced  data,  s2  is  given  by 


(71) 
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0 


2  AS (45)  +  AS  (135) 

In  this  case,  the  numerator  for  both  the  undifferenced  and  differenced  data  is 
approximately  zero,  but  the  denominator  is  much  smaller  for  the  differenced  data.  So  to 
any  extent  that  the  numerator  is  not  zero,  the  differenced  calculation  of  s2  will  always  be 

larger.  A  non-zero  numerator  can  result  from  noise  in  the  measurements  and  polarizer 
misalignment. 
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Figure  29.  Images  showing  the  effect  on  the  Stokes  elements  and  DoLP  of  differencing  the  high  and 
low  temperature  data  sets.  The  external  and  internal  polarizers  are  aligned  horizontally  which 
should  result  in  a  value  of  sx  close  to  one  and  s2  close  to  zero  with  a  DoLP  close  to  100%.  In  the 
differenced  data  this  is  true  whereas  in  the  undifferenced  data  this  is  far  from  true  due  to  the 

reflected  and  emitted  components. 
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Analysis  of  Results 

Evaluation  of  System  Matrices 

Comparing  the  experimentally  determined  system  matrix  to  the  ideal  system 
matrix  provides  valuable  insight  into  how  ideal  or  non-ideal  the  system  is  performing. 

All  of  the  non-ideal  effects  can  be  captured  in  a  single  Mueller  matrix  as  follows: 

Lm  =  W  •  5  =  WidealMA  •  S  (73) 

where  the  measured  system  matrix  has  been  split  into  the  ideal  system  matrix  multiplied 
by  a  Mueller  matrix,  MA ,  called  the  Mueller  deviation  matrix  [20].  From  Equation  (73), 

Ma  can  be  determined  with 

Ma=W,;eaW.  (74) 

If  the  system  were  ideal,  this  matrix  would  be  equal  to  the  identity  matrix,  so  any 
deviations  from  the  identity  matrix  is  an  indication  of  non-ideal  effects.  These  non-ideal 
effects  are  a  result  of  the  non-ideal  nature  of  the  instrument  polarizer,  the  non-zero 
diattenuation  of  the  HSI  instrument  following  the  polarizer,  the  non-linearity  of  the 
detector,  and  other  effects  such  as  noise.  The  Mueller  deviation  matrix  also  has  the 
interesting  property  that  it  converts  an  input  Stokes  vector  into  the  Stokes  vector  that 
would  be  measured  by  the  system  if  it  were  assumed  to  be  ideal,  i.e. 

Sm=MA-Sm.  (75) 

This  provides  another  useful  way  to  understand  physically  what  the  system  is  doing.  By 
multiplying  the  Mueller  deviation  matrix  by  various  ideal  input  Stokes  vectors  and  seeing 
what  comes  out,  one  gains  an  understanding  of  how  well  the  system  is  performing  and 
what  the  different  biases  are.  In  fact,  another  method  of  calibrating  the  instrument  via 
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direct  determination  of  the  Mueller  deviation  matrix  can  be  developed.  Solving  for  the 
Mueller  deviation  matrix  with  multiple  input  Stokes  vectors  gives 

MA=Sm-S;.  (76) 

Once  the  Mueller  deviation  matrix  is  determined,  Equation  (73)  can  be  used  to  determine 
the  system  matrix  via 

W  =  W]dealMA  (77) 

and  the  data  reduction  matrix  can  then  be  determined  by  taking  the  pseudo-inverse  of  W . 
This  method  is  effectively  the  same  as  the  previous  method  mathematically  and  is  simply 
provided  as  an  additional  option. 

Recall  that  the  ideal  system  matrix  is  given  by 

'11  0^ 

WkM=  !  q  (78) 

(l  0  -lj 

and  as  mentioned  above,  the  corresponding  ideal  Mueller  deviation  matrix  is  given  by  the 
identity  matrix,  i.e. 

'1  0  0^ 

MA  =  0  1  0  .  (79) 

1°  0  lj 

Using  the  experimental  setup  and  methodology  described  above,  a  system  matrix  was 
computed  for  the  Telops  Hypercam.  Since  the  Hypercam  is  an  imaging  spectrometer, 
this  had  to  be  done  on  a  pixel-by-pixel  and  band-by-band  basis,  creating  a  data  cube  of 
system  matrices.  Since  it  is  impractical  to  report  the  results  for  each  pixel  at  each  band,  a 
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spectrally  and  spatially  averaged  system  matrix  was  also  computed.  The  resulting  system 


matrix  is 


1.015 

1.015 

0.010 

1.000 

0.011 

0.999 

0.992 

-0.988 

0.018 

0.992 

0.010 

-1.004 

(80) 


While  not  perfectly  matching  the  ideal  matrix,  the  results  show  that  on  an  averaged  basis, 
the  system  is  performing  very  close  to  ideally.  The  fact  that  the  first  column  is  not  a 
column  of  ones  suggests  that  the  radiometric  calibration  was  not  perfect  or  that  the 
temperature  of  the  blackbody  sources  were  not  stable.  To  prevent  applying  a  polarimetric 
bias  to  unpolarized  light,  the  first  row  must  be  set  to  ones.  The  corresponding  average 
Mueller  deviation  matrix  is 


M 


A 


(1.000  0.007  0.001  N 
0.000  1.001  -0.014 
v0.000  0.011  1.001  , 


(81) 


As  mentioned  above,  the  Mueller  deviation  matrix  is  used  to  calculate  the  Stokes  vector 
that  the  system  would  measure  without  any  correction  and  this  provides  another  good 
way  to  evaluate  how  well  the  system  is  performing  by  multiplying  this  Mueller  deviation 
matrix  by  various  ideal  Stokes  vectors.  For  instance,  if  perfectly  horizontally  polarized 
light  were  incident  on  the  instrument,  the  output  Stokes  measurement  would  be 


= 


'1.000  0.007  0.00P 

(l) 

"1.007) 

0.000  1.001  -0.014 

1 

= 

1.013 

v0.000  0.011  1.001  , 

w 

v0.015y 

(82) 


Similarly,  vertically  polarized  incident  light  gives 
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These  results  indicate  less  than  1.5%  error  in  most  cases  before  any  further  polarimetric 
calibration  is  performed.  The  DoLPs  associated  with  these  four  measurements  are 
99.4%,  100.9%,  100.1%,  and  100.2%  respectively.  While  DoLPs  of  over  100%  are 
unphysical,  recall  the  mathematical  argument  in  Equation  (68)  and  the  fact  that  the 
system  matrix  is  computed  with  a  least  squares  solution  to  measurements  subject  to  noise 
and  limited  input  accuracy. 

Since  the  instrument  is  an  imaging  spectrometer,  it  is  also  important  to  evaluate 
the  performance  on  a  pixel-by-pixel  and  band  by  band  basis.  To  start  with,  four  random 
pixels  at  four  random  bands  were  selected  to  report  the  system  matrix  of  those  individual 
pixels  at  the  respective  individual  bands.  The  results  are  given  below: 
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^1.027 

1.023 

-0.010^ 

"1.035  1.025 

-0.012^ 

w,  = 

1.004 

0.005 

1.010 

W,  = 

1.005  -0.002 

1.001 

0.991 

-0.987 

0.016 

T  T  2 

0.968  -0.989 

0.014 

,0.979 

-0.001 

-0.998 j 

0.992  -0.026 

-1.022, 

(46,44,1049) 

(48,54,998) 

"1.014 

1.036 

0.004 

"1.032  1.041 

-0.008^ 

W3  = 

1.000 

0.007 

0.967 

w4  = 

1.003  0.021 

1.003 

0.990 

-1.000 

0.031 

0.990  -0.972 

0.020 

v0.996 

0.014 

-1.021, 

v0.985  -0.012 

-0.981, 

(49,51,1101)  (54,40,915) 


where  the  numbers  in  parenthesis  below  the  matrices  corresponds  to  the  pixels’  row, 
column,  and  band  in  wavenumbers,  respectively.  This  shows  that  the  system  is  generally 
within  about  4%  of  the  ideal  case  for  the  pixels  selected  even  at  the  individual  pixel  and 
individual  band  level  when  no  spatial  or  spectral  averaging  has  taken  place.  A  more 
complete  view  of  the  spectral  nature  of  each  of  the  four  pixels  selected  is  shown  in  Figure 
30.  The  elements  of  the  system  matrix  for  each  of  the  four  pixels  discussed  above  are 
plotted  spectrally.  While  there  is  some  noise  in  the  spectral  elements,  the  results  are  good 
overall,  showing  that  the  elements  of  the  system  matrix  are  spectrally  flat  for  each 
individual  pixel.  The  lines  are  color  coded  as  a  visual  aid.  Blue  lines  correspond  to 
elements  with  a  value  near  one,  black  lines  are  near  zero,  and  red  lines  are  near  negative 
one.  It  is  common  to  do  some  spatial  averaging  over  a  group  of  pixels  that  are  imaging  a 
spatially  uniform  part  of  the  scene.  Any  averaging  of  this  nature  will  reduce  the  noise 
and  flatten  the  lines  out  further.  Figure  31  shows  the  plot  of  the  spatially  averaged 
system  matrix  with  the  standard  deviation  represented  with  shaded  bands.  Note  that  the 
scale  on  this  plot  is  much  narrower  than  the  scale  on  the  individual  pixel  plots  to  aid 
visualization  of  the  noise  in  the  measurements.  It  is  clear  that  the  majority  of  the  noise  is 
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removed  as  a  result  of  the  spatial  averaging.  The  average  standard  deviation  in  the 


measurement  of  the  system  matrix  ranges  from  0.35%  to  0.77%. 
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a.)  Spectral  plot  of  the  W  matrix  elements  for  pixel  1. 
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b.)  Spectral  plot  of  the  W  matrix  elements  for  pixel  2. 
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c.)  Spectral  plot  of  the  W  matrix  elements  for  pixel  3. 
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d.)  Spectral  plot  of  the  W  matrix  elements  for  pixel  4 
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Figure  30.  Spectral  plots  of  the  system  matrix  elements  for  each  of  the  four  pixels  discussed  above. 
The  elements  of  W  for  each  individual  pixel  are  spectrally  flat  other  than  the  random  noise  seen  in 
the  plots.  The  lines  are  color  coded  as  a  visual  aid.  Blue  lines  correspond  to  elements  with  a  value 
near  one,  black  lines  are  near  zero,  and  red  lines  are  near  negative  one. 
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Figure  31.  Spectral  plot  of  the  system  matrix,  W,  after  spatial  averaging.  The  standard  deviation  is 
represented  with  the  shaded  bands  around  measured  values.  Since  this  is  a  spatial  average,  the  noise 
is  reduced  considerably.  The  scale  on  this  plot  is  zoomed  in  to  allow  better  visualization  of  the  noise 

in  the  measurements. 


Another  interesting  result  is  shown  in  Figure  32  where  each  element  of  the  band- 
averaged  Mueller  deviation  matrix  is  plotted  for  each  pixel  creating  a  spatial  map.  The 
diagonal  elements  are  all  close  to  one  and  the  off-diagonal  elements  are  all  close  to  zero 
as  expected.  The  checkerboard  pattern  seen  in  the  off-diagonal  elements  is  believed  to  be 
due  to  non-linear  effects  in  the  read  out  integrated  circuitry  (ROIC)  and  represents  a 
current  limit  on  data  accuracy.  Calibrating  this  effect  out  is  the  subject  of  a  separate, 
ongoing  research  effort. 
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Figure  32.  Each  element  of  the  band-averaged  Mueller  deviation  matrix  is  plotted  for  each  pixel. 
The  diagonal  elements  are  all  close  to  one  and  the  off-diagonal  elements  are  all  close  to  zero  as 
expected.  The  checkerboard  pattern  in  the  off-diagonal  elements  is  believed  to  be  due  to  non-linear 
effects  in  the  read  out  integrated  circuitry  (ROIC)  and  represents  a  current  limit  on  data  accuracy. 


Calibration  Verification  on  an  Unpolarized  Source 

Now  that  the  data  reduction  matrix  has  been  determined  experimentally,  it  is 

necessary  to  verify  how  well  the  calibration  is  working  under  various  input  states.  From 

here  on,  the  term  “uncorrected  data”  refers  to  data  that  has  been  processed  and 

polarimetrically  calibrated  under  the  assumption  of  an  ideal  polarizer,  and  the  term 

“corrected  data”  refers  to  the  additional  correction  for  the  non-ideal  polarizer  by 

application  of  the  experimentally  determined  data  reduction  matrix.  The  first  verification 

will  measure  the  system  response  to  unpolarized  light.  It  is  important  to  note  here  that 
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the  diattenuation  of  the  polarizers  has  no  effect  on  the  system  response  to  unpolarized 
light.  This  can  be  verified  by  plugging  Sj  =  0 ,  s2  -  0 ,  and  s3  =  0  into  Equation  (45)  (See 

Chapter  V  for  more  details  on  this  equation  for  the  non-ideal  polarizer  case).  The 
correction,  therefore,  should  have  no  effect  on  the  measured  Stokes  vector.  To 
demonstrate  this,  a  blackbody  source  was  used  as  an  unpolarized  source  to  determine 
whether  or  not  the  system  showed  any  polarimetric  bias.  A  6-inch  blackbody  was  used  to 
flood  fill  the  detector  and  200  data  cubes  were  collected  for  each  instrument  polarizer 
angle.  On-board  blackbody  calibration  data  was  also  collected  at  each  polarization  angle. 
Each  of  the  200  data  cubes  at  each  angle  was  then  individually  calibrated  and  used  to 
determine  the  Stokes  vector  of  the  incoming  light.  This  was  done  spectrally  on  a  pixel- 
by-pixel  basis.  Then,  the  spatial  and  band  average  of  each  data  cube  was  taken.  This 
generated  an  averaged  Stokes  vector  for  each  of  the  200  measurements.  Next,  the 
histogram  of  the  200  Stokes  vectors  was  determined  and  the  normalized  sx  and  s2 
distributions  were  plotted.  For  each  histogram,  a  normal  distribution  curve  was  fitted  to 
the  data.  The  results  are  shown  in  Figure  33.  The  plot  shows  that  the  mean  of  s{  and 

s2  are  8.6 x  10  4  and  7.6x10  4,  respectively,  for  the  uncorrected  data.  The  corrected  data 
is  also  plotted  and  lines  up  almost  on  top  of  the  uncorrected  data  as  expected.  Ideally, 
both  sx  and  s2m  would  be  zero  for  unpolarized  light.  In  this  case,  there  appears  to  be  a 

very  small  positive  bias  on  both  slm  and  s2m  .  It  is  not  clear  exactly  what  is  causing  the 

bias;  however,  some  possibilities  include  non-stable  target  or  calibration  blackbody 

sources  or  the  non-linearity  of  the  detector.  As  mentioned  before,  the  6-inch  blackbody 

source  has  a  temperature  stability  of  ±  0.5  °C.  To  determine  if  this  was  a  possible 
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candidate  for  the  bias,  a  simulation  was  run  to  determine  the  temperature  drift  that  would 
be  required  between  measurements  at  orthogonal  polarizer  angles  to  produce  this  bias. 
The  results  of  the  simulation  showed  that  a  temperature  difference  of  0.5  °C  between 
orthogonal  measurements  produces  a  band-averaged  bias  of  approximately  8.6  xlCT4 . 
This  is,  therefore,  a  feasible  source  of  the  bias.  Since  this  is  an  important  source  of 
uncertainty  introduced  by  the  experimental  setup,  the  thermal  stability  of  the  6-inch 
blackbody  was  independently  verified  rather  than  relying  on  the  manufacturer’s 
specification.  The  results  in  Figure  34  show  that  the  6-inch  blackbody  is  indeed 
oscillating  approximately  0.5  °C  about  the  set  point  every  1.2  minutes. 
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Figure  33.  Histogram  of  Si  and  s2  Stokes  elements  for  200  measurements  of  blackbody  radiation. 
The  histograms  are  fitted  with  a  normal  distribution  curve.  In  the  ideal  case,  the  mean  of  Si  and  s2 
would  be  zero.  This  shows  that  there  is  a  small  bias  in  the  measurements,  possibly  to  unstable 
blackbody  sources  or  the  non-linearity  of  the  detector. 
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Figure  34.  Temperature  stability  of  the  6-inch  blackbody  over  a  period  of  just  more  than  6  minutes. 
This  verifies  that  the  temperature  is  oscillating  approximately  0.5  °C  with  a  period  of  less  than  1.5 

minutes. 


Calibration  Verification  on  Polarized  Sources 

Now,  how  well  this  experimentally  determined  matrix  is  able  to  correct  for  the 
non-idealities  of  the  system  is  explored  when  polarized  light  is  input.  The  band-averaged 
extinction  ratio  of  the  Hypercam  polarizer  is  approximately  400.  This  gives  a 
diattenuation  of 


D  ,  =  1  =  400  1  =  0.995  . 


pol 


R  + 1  400  +  1 


This  means  that  the  uncorrected  polarimetric  calibration  should  achieve  results  on  the 
order  of  0.5%  from  the  ideal  case  if  noise  and  other  uncertainties  have  been  sufficiently 
reduced.  To  evaluate  the  effectiveness  of  the  correction  on  polarized  light,  the  measured 
Stokes  parameters  output  from  the  system  were  compared  to  the  input  Stokes  vectors  for 
both  the  uncorrected  case  and  the  corrected  case.  The  spatially  and  spectrally  averaged 


80 


results  are  shown  in  Figure  35.  In  each  plot,  the  solid  line  shows  the  input  Stokes  value 
as  a  function  of  external  polarizer  orientation  and  the  shaded  band  represents  the 
uncertainty  in  the  external  polarizer  orientation.  The  blue  markers  are  for  the  uncorrected 
data  and  the  red  markers  are  for  the  corrected  data.  The  ‘x’  markers  are  at  the  calibration 
points  and  the  ‘+’  markers  are  in-between  calibration  points.  It  is  difficult  to  visually  tell 
the  difference  between  the  corrected  and  uncorrected  data  from  the  plots,  which  again 
confirms  that  the  polarimetric  calibration  under  the  ideal  polarizer  assumption  works  very 
well.  To  see  the  effect  of  the  correction,  a  zoomed  in  view  is  shown  for  a  specific  data 
point.  Now  it  can  be  seen  that  the  corrected  data  is  closer  to  the  input  line  than  the 
uncorrected  data  for  that  specific  data  point. 

To  visualize  the  spectral  behavior  of  the  measured  Stokes  values,  both  on  a 
spatially  averaged  and  a  pixel-by-pixel  basis,  Figure  36  shows  a  3-D  plot  of  the  data 
matrix.  The  two  plots  on  the  left  show  the  measured  values  of  .v,  and  s2  for  one  pixel 

and  the  plots  on  the  right  show  the  measured  values  of  s1  and  s2  for  a  spatial  average  of 

317  pixels.  The  results  indicate  that,  while  there  is  some  noise  in  the  measured  values  of 
one  pixel,  the  measured  Stokes  values  are  flat  spectrally.  This  is  clear  on  the  averaged 
data.  Note  that  Figure  36  only  shows  the  uncorrected  data;  however,  as  in  Figure  35,  it  is 
hard  to  discern  any  visual  difference  between  the  corrected  and  uncorrected  data  in  these 
plots. 
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Figure  35.  Comparison  of  the  measures  Stokes  values  for  the  corrected  and  uncorrected  data  for  a.) 
Si  and  b.)  s2.  The  solid  line  is  the  input  Stokes  values  with  a  shaded  band  representing  the 
uncertainty  in  the  external  polarizer  orientation. 
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Figure  36.  Shows  the  spectral  nature  of  each  of  the  12  measurements  for  a  single  pixel  (left)  and  for 
317  spatially  averaged  pixels  (right).  The  top  plots  are  Si  and  the  bottom  plots  are  s2. 


To  further  analyze  the  performance  of  the  correction,  the  error  between  the 
measured  data  and  the  input  values  are  calculated  for  both  the  corrected  and  uncorrected 
data.  In  this  case,  the  data  at  the  calibration  points  are  plotted  separately  from  the  data 
taken  in-between  calibration  points.  This  helps  to  understand  how  much  error  is  in  the 
least  squares  solution  at  the  calibration  points.  Figure  37  shows  these  results.  In  all 
cases,  the  correction  reduced  the  average  error  although  there  are  still  points  where  the 
corrected  data  has  a  greater  error  than  the  uncorrected  data.  The  data  shows  that  on  a 
spatially  and  band  averaged  basis,  the  correction  improves  the  error  by  0.5%  or  less  and 
that  there  is  a  residual  error  ranging  from  approximately  1-2%  even  after  calibration.  The 
primary  cause  of  the  residual  error  is  likely  to  be  the  thermal  drift  of  the  6-inch  blackbody 
as  discussed  in  the  previous  section.  The  correction  of  0.5%  is  consistent  with  the  fact 
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that  the  instrument  polarizer  is  within  about  0.5%  of  an  ideal  polarizer  based  on  its 
diattenuation. 
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Figure  37.  a.)  Error  between  the  measured  data  and  the  input  data  both  before  and  after  correction 
at  the  12  calibration  points,  b.)  Error  for  the  data  collected  in-between  the  calibration  points. 

Next,  the  same  error  plots  are  shown  in  Figure  38  for  an  individual  pixel  at  two 
individual  bands,  one  near  the  middle  of  the  bandwidth  and  one  near  the  edge  of  the 
bandwidth  where  the  SNR  is  slightly  decreased  due  to  the  sensor's  sensitivity.  As 
expected,  the  average  error  increases,  now  ranging  from  roughly  1-4%;  however,  on 
average,  the  error  after  correction  is  less  than  the  error  before  correction  for  both  s1  and 

s2 .  To  gain  a  better  understanding  of  the  error  of  one  pixel  across  the  entire  bandwidth. 
Figure  38  shows  the  error  of  the  same  pixel  plotted  spectrally  for  each  of  the  12  measured 
Stokes  vectors  before  and  after  correction  for  .v,  and  s2 .  Additionally,  since  there  will 

likely  be  some  spatial  averaging  done  during  post-processing  for  a  typical  scene,  Figure 
39  shows  the  spectral  error  for  317  spatially  averaged  pixels.  The  results  show  that  on  an 
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individual  pixel  basis,  the  error  of  the  corrected  data  tends  to  stay  either  above  the 
uncorrected  error  for  the  entire  band  or  below  the  uncorrected  error  for  the  entire  band 
without  jumping  back  and  forth  too  much.  On  a  spatially  averaged  basis  this  trend  is 
even  stronger  since  the  error  flattens  out  considerably  across  the  band.  The  average  error 
across  the  band  is  also  reduced.  In  order  to  aid  visual  comparison  between  the  individual 
pixel  results  and  the  averaged  pixel  results,  they  were  both  plotted  on  the  same  scale. 
Notice  that  in  all  of  the  error  plots  shown  below,  there  appears  to  be  a  sinusoidal 
oscillatory  nature  to  error  magnitude.  It  is  speculated  that  this  oscillatory  nature  in  the 
error,  is  again  related  to  the  oscillatory  thermal  fluctuations  of  the  6-inch  blackbody 


shown  in  Figure  34. 


a.) 
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Figure  38.  Error  of  the  measured  data  for  an  individual  pixel  at  a.)  near  the  center  of  the  bandwidth 
and  b.)  near  the  edge  of  the  bandwidth  where  detector  sensitivity  is  slightly  decreased. 
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Figure  39.  Shows  the  measured  spectral  error  for  one  pixel  for  a.)  sx  and  b.)  s2  for  each  of  the  12 
input  Stokes  vectors.  The  label  to  the  right  of  each  plot  indicates  the  orientation  of  the  external 

polarizer. 
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Figure  40.  Shows  the  measured  spectral  error  of  317  spatially  averaged  pixels  for  a.)  Si  and  b.)  s2  for 
each  of  the  12  input  Stokes  vectors.  The  label  to  the  right  of  each  plot  indicates  the  orientation  of  the 

external  polarizer. 

Noise  and  Error  Analysis 

Next,  it  is  important  to  quantify  the  noise  in  the  measurements  in  order  to  gain 
further  insight  into  the  error  results  shown  above.  Theoretically,  if  the  noise  in  the 
measurements  was  zero  and  the  detector  was  perfectly  linear,  then  the  polarimetric 
calibration  correction  should  perfectly  correct  the  measured  Stokes  vectors  back  to  the 
known  input  Stokes  vectors.  The  more  random  noise  in  the  measurements  and  the  more 
non-systematic  errors,  the  worse  the  least  squares  solution  is  going  to  be  and  the  worse 
the  polarimetric  calibration  is  going  to  perform.  Unfortunately,  the  SNR  for  this 
experiment  is  forced  down  by  the  fact  that  the  hot  and  cold  collections  are  differenced. 
This  differencing  effectively  reduces  the  signal  while  increasing  the  random  noise. 
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Some  typical  noise  sources  include  photon  shot  noise,  shot  noise  from  the  dark 
current,  Johnson  noise,  preamplifier  noise,  and  quantization  noise  [28].  Additional 
sources  of  noise  that  are  introduced  by  the  limitations  of  the  experimental  setup  include 
random  error  when  setting  the  external  polarizer  angle,  random  fluctuations  of  the 
external  blackbody  sources,  and  random  or  inconsistent  thermal  fluctuations  in  the 
external  polarizer  between  hot  and  cold  measurements.  In  addition  to  these  random 
noises,  systematic  errors  can  exist  which  will  bias  the  data  in  the  same  way  for  each 
measurement.  Sources  of  possible  systematic  error  include  the  uncertainty  in 
determining  the  absolute  zero-angle  of  the  external  polarizer  and  the  uncertainty  in  the 
spectral  diattenuation  of  the  external  polarizer.  Thermal  drift  of  the  external  polarizer 
could  also  be  considered  a  systematic  error  if  the  thermal  drift  was  constant  and  affected 
each  measurement  the  same  way.  A  systematic  error  in  the  external  polarizer  angle  will 
result  in  a  horizontal  shift  of  the  input  Stokes  vector  in  Figure  35  and  other  similar  plots, 
and  will  cause  the  calculated  error  before  correction  to  be  artificially  high.  An  error  in 
the  diattenuation  of  the  external  polarizer  will  result  in  a  vertical  scaling  of  the  input 
Stokes  vector  and  will  cause  the  corrected  Stokes  values  to  be  scaled  by  the  ratio  of  the 
measured  diattenuation  to  the  true  diattenuation.  Thermal  drift  of  the  external  polarizer 
between  hot  and  cold  measurements  could  cause  artificially  high  or  low  measured  values 
for  s1  and  s2 .  For  the  error  analysis,  the  relative  systematic  error  due  to  an  incorrect 

determination  of  the  external  polarizer  zero-angle  is  assumed  to  be  zero  since  it  was 
removed  by  fitting  a  sinusoid  to  the  uncorrected  measurements  and  correcting  for  the 
phase  shift  between  the  input  Stokes  vector  and  the  fitted  sinusoid.  Therefore,  the 
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relative  error  before  the  correction  is  a  result  of  random  error/noise  and  the  assumption  of 
an  ideal  polarizer.  The  correction  should  then  remove  the  assumption  of  an  ideal 
polarizer  and  leave  only  error  due  to  random  noise,  detector  non-linearity,  and  the 
thermal  fluctuations  of  the  6-in  blackbody  discussed  before.  This  means  that  the  random 
noise  in  the  measurements  of  the  input  Stokes  vectors  (after  frame  averaging)  must  be 
smaller  than  the  error  due  to  the  assumption  of  an  ideal  polarizer  in  order  to  have  a 
definitive  effect.  If  this  is  not  true,  then  the  uncertainty  in  the  solution  of  the  system 
matrix  will  be  larger  than  the  0.5%  correction  it  is  attempting  to  make.  It  was  already 
shown  in  Figure  3 1  that  the  noise  in  the  measurement  of  the  system  matrix  elements 
ranged  from  0.35-0.77%.  Therefore,  it  will  be  difficult  to  make  a  0.5%  correction  unless 
the  noise  is  further  reduced.  Additionally,  the  noise  statistics  of  s1  and  s2  must  be  used 

to  compare  to  the  error  reported  for  s1  and  s2  earlier.  This  is  shown  in  Figure  41.  The 

noise  ranges  from  around  1-2%.  This  noise  can  account  for  some  of  the  residual  error 
after  calibration. 
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Figure  41.  Spectrally  resolved  noise  of  the  Si  and  s2  measurements  when  viewing  horizontally 
polarized  light.  The  noise  was  determined  by  taking  the  standard  deviation  of  the  mean  for  the  64 

frames  acquired  at  each  band. 


NESR  and  NESDoLP  of  an  Unpolarized  Source 

In  addition  to  trying  to  correct  for  the  non-ideal  nature  of  the  instrument  polarizer, 
there  was  also  a  desire  to  understand  what  the  noise  equivalent  spectral  radiance  (NESR) 
and  noise  equivalent  DoLP  (NEDoLP)  of  the  system  are. 

To  determine  the  NESR  from  experimental  data,  a  45  °C  blackbody  was  placed  in 

front  of  the  Telops  so  that  the  full  detector  array  was  flood  filled.  The  integration  time 

was  set  to  200  jus  and  the  spectral  resolution  was  16  cm"1.  The  blackbody  source 

provides  the  ideal  condition  to  measure  the  NESR  since  the  temporal  variations  of  the 

scene  radiance  are  negligible  with  a  well-controlled  blackbody  [3].  Next,  200  data  cubes 

were  collected  at  each  of  the  four  standard  instrument  polarizer  angles  (0°,45o,90°,  and 

135°).  A  radiometric  calibration  was  performed  at  each  of  the  instrument  polarizer  angles 

and  bad  pixels  were  removed.  The  NESR  is  then  determined  for  each  pixel  by 
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calculating  the  temporal  standard  deviation  of  the  measured  radiance  at  each  waveband. 
The  pixel  average  of  the  measured  NESR  for  each  instrument  polarizer  angle  is  shown  in 
Figure  43.  The  NESR  values  measured  are  consistent,  although  slightly  lower,  than  the 
NESR  values  reported  by  Telops  in  [3].  It  makes  sense  that  the  values  are  lower  since  the 
values  reported  by  Telops  were  for  a  4  cm'1  resolution  as  opposed  to  a  16  cm'1  resolution 
and  NESR  increases  with  higher  spectral  resolution  [2],  It  is  believed  that  the  difference 
between  the  NESR  at  0°  and  90°  is  due  to  the  non-linearity  of  the  detector  response.  This 
difference,  however,  had  little  effect  on  the  calibration’s  ability  to  normalize  the  system 
response  since  the  SNR  is  greater  than  100  across  the  band.  This  can  be  shown  in  Figure 
42  where  the  radiance  measured  by  the  detector  at  each  instrument  polarizer  angle  is 
plotted.  There  is  very  little  deviation  between  the  measured  radiance  at  each  polarizer 
angle. 


Figure  42.  Measured  NESR  of  the  Telops  instrument  at  each  polarizer  angle. 
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Figure  43.  Measured  radiance  when  viewing  an  unpolarized  blackbody  source  at  each  instrument 

polarizer  angle  after  calibration. 


Next,  the  NESDoLP  is  determined.  The  NESDoLP  is  effectively  the  minimum 
DoLP  the  system  can  discern  and  is  based  on  the  noise  in  the  measurement  of  5,1  and  s2 . 
Recall  that  DoLP  is  given  by 


DoLP  =  y[sf+sf ,  (86) 

where  in  this  case,  s1  and  s2  have  already  been  normalized  by  ,v0 .  This  means  that  the 
accuracy  in  the  measurement  of  DoLP  can  never  be  better  than  the  propagated 
uncertainty  in  st  and  s2 .  The  propagation  of  uncertainty  of  Equation  (86) ,  and 
therefore  the  NEDoLP,  is  given  by 


NEDoLP  = 


\(  oDoLP 


v 


eJSj 

J 


+ 


dDoLP 

ds~ 


y 

Ss2 

J 


(87) 


which  gives 


94 


NEDoLP= 

where  Sst  and  Ss2  are  the  temporal  standard  deviation  of  s1  and  s2 .  Doing  this  at  each 

pixel  and  at  each  waveband  generates  the  spectrally  resolved  NEDoLP,  or  the  noise 
equivalent  spectral  DoLP  (NESDoLP)  given  by 

NESDoLP  (i,j,cr)= 

where  i  is  the  pixel  row  index,  j  is  the  pixel  column  index,  and  cr  is  the  spectral 


(/,;,<t)  •£.?,(/,  7,  cr) 
Js;(i,j,a)  +  s2(i,j,a) 


s2(i,j,c>)-8s2(i,j,c>) 
Jsf  (i,  j,  cr)  +  .S' 2  (i, j,  cr) 


dependence  on  wavenumbers. 

Further  insight  into  the  actual  system  parameters  that  determine  NEDoLP  is 
provided  by  Jones  and  Persons.  In  their  paper  (Ref  [27]),  they  derive  a  model  for  the 
NEDoLP  of  a  micro-grid  polarimeter.  This  model  is  useful  in  understanding  the  system 
parameters  which  affect  the  NEDoLP.  Their  model  for  NEDoLP  is  given  by 


where  8  is  the  ratio  of  read  noise  to  shot  noise,  /?  is  the  ratio  of  signal  electrons 


accumulated  in  the  well  to  total  electrons  accumulated  in  the  well,  D  is  the  diattenuation 
of  the  polarizer,  and  n  is  the  well  capacity.  From  this,  the  following  observations  can  be 
made:  First,  if  the  read  noise  is  significant  compared  to  the  shot  noise,  then  NEDoLP  is 
degraded.  Second,  the  NEDoLP  is  improved  with  a  larger  well  size  of  the  FPA  so  the 
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polarimetric  performance  will  benefit  from  a  large  well  capacity.  Third,  a  better  SNR 
will  improve  NEDoLP.  Finally,  a  more  ideal  polarizer  will  improve  the  NEDoLP. 

In  order  to  experimentally  determine  the  NESDoLP,  the  same  dataset  from  the 
NESR  determination  was  used.  The  NESDoLP  was  calculated  for  each  pixel  at  each 
waveband  using  Equation  (89).  This  was  done  for  both  the  corrected  and  uncorrected 
data.  As  demonstrated  in  the  "Calibration  Verification  on  Unpolarized  Source"  section,  if 
the  detector  non-linearity  is  small  and  the  calibration  sources  are  stable,  the  correction 
should  have  little  effect  on  unpolarized  light,  and  therefore  the  NESDoLP  should  be  the 
same  for  the  corrected  and  uncorrected  data.  The  pixel  average  of  the  NESDoLP  for  the 
corrected  and  uncorrected  data  is  shown  in  Figure  44.  As  expected,  the  corrected  and 
uncorrected  data  give  the  same  results.  The  NESDoLP  ranges  from  0.32%  to  0.55% 
across  the  band.  There  are  a  few  things  to  keep  in  mind  here.  First,  this  is  the  NESDoLP 
for  one  measurement.  If  N  measurements  are  averaged  then  the  NESDoLP  is  improved 
by 


NESDoLP'= 


NESDoLP 

s[N 


(91) 


Second,  the  NESDoLP  will  change  depending  on  the  SNR  in  the  measurements.  Third, 
this  assumes  that  all  other  error  has  been  removed  and  that  the  only  remaining  uncertainty 
in  the  measurement  comes  from  the  random  noise.  Assuming  that  the  blackbody  is 
randomly  polarized,  the  calculated  DoLP  should  be  approximately  the  same  as  the 
NESDoLP.  If  this  is  not  true,  then  some  other  error  or  uncertainty  other  than  random 
noise  remains  in  the  measurements.  Figure  45  shows  the  spectral  DoLP  plotted  with  the 
NESDoLP.  Since  the  measured  DoLP  is  larger  than  the  NESDoLP,  there  appears  to  be 
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some  uncertainty  in  the  measurements  which  has  not  been  removed.  The  figure  shows 
that  on  average  the  DoLP  errors  are  approximately  0.12%  larger  than  what  the  NESDoLP 
predicts.  This  again,  could  be  due  to  blackbody  sources  that  are  fluctuating  slightly  or  it 
could  be  due  to  the  non-linearity  of  the  detector.  It  has  been  shown  that  the  ability  to 
measure  low  DoLP  signatures  is  greatly  influenced  by  the  accuracy  of  the  2-point 
calibration  [29],  [30]  [31],  [25], 
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Figure  45.  DoLP  and  NESDoLP  plotted  together.  The  DoLP  error  is  approximately  0.12%  larger 

than  predicted  by  the  NESDoLP. 


NESR  and  NESDoLP  of  a  Polarized  Source 

Now,  the  exact  same  plots  as  above  are  plotted  for  a  source  that  is  close  to 

horizontally  polarized  (after  fitting  a  sinusoid  to  the  data,  it  was  actually  determined  to  be 

polarized  at  1.57°).  The  horizontally  polarized  source  was  generated  using  the  setup 

shown  in  Figure  22  and  the  same  differencing  method  described  earlier.  One  expectation 

here  would  be  that  the  correction  should  now  have  an  effect  since  polarized  light  is  input 

into  the  system.  Another  expectation  would  be  that  the  noise  has  more  of  an  effect  on  the 

data  because  the  SNR  is  much  lower  due  to  the  differencing.  Figure  46  shows  the  NESR 

at  each  instrument  polarizer  angle.  The  NESR  is  roughly  the  same  as  it  was  for  the 

unpolarized  source.  Figure  47  shows  the  measured  radiance  at  each  instrument  polarizer 

angle.  Notice  that  because  of  the  differencing,  the  radiance  measured  at  each  angle  is 

small  compared  to  the  radiance  measured  for  the  unpolarized  source.  Also,  the  radiance 
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at  45°  and  135°  are  off  slightly  due  to  the  1.57°  offset  mentioned  before.  Figure  48  shows 
the  measured  DoLP  before  and  after  correction.  This  time  the  correction  improves  the 
DoLP  measurement  by  an  average  of  0.56%. 


Figure  46.  Measured  NESR  for  horizontally  polarized  light  at  each  instrument  polarizer  angle. 
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DoLP  [%]  Radiance  [|iW/cm2sr  cm 


Figure  47.  Measured  radiance  at  each  instrument  polarizer  angle  for  a  horizontally  polarized 

source. 


Figure  48.  Measured  DoLP  before  and  after  correction  for  horizontally  polarized  light.  This  time 
the  correction  improves  the  measurement  by  0.56%  on  average. 
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V.  Polarimetric  Performance  Modeling 


This  chapter  utilizes  the  mathematical  frame  work  developed  in  the  previous 
chapter  to  analyze  the  effect  of  having  a  non-ideal  instrument  polarizer  and  a  non-ideal 
HSI  system  Mueller  matrix.  In  the  previous  chapter,  the  Mueller  matrix  of  an  ideal 
polarizer  was  used  to  determine  the  functional  form  of  m01  (0  )  and  m02  (0  )  in  Equation 


(44).  The  functional  form  of  these  elements  must  now  be  determined  using  the  more 
general  Mueller  matrix  of  a  linear  polarizer  that  includes  diattenuation.  This  Mueller 


matrix  is  given  by 


2 


q  +  r 
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(q-r)  sin  20 p 
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2  yfqr) 


where  q  is  the  maximum  transmission  of  linearly  polarized  light  and  r  is  the  minimum 


transmission  of  linearly  polarized  light.  Diattenuation  is  then  given  by  (q-r)/(q  +  r ) 
and  the  extinction  ratio  is  given  by  ql  r .  Plugging  this  into  Equation  (41)  and  noting  that 
sin  20  cos  20  is  always  zero  at  the  angles  used  in  the  calibration  (i.e.  0°,  45°,  90°,  and 
135°),  it  is  found  that 


m(n  (0„ )  = 


(q  -  r ) cos 26 p  +  m'0l  i^(q  +  r ) cos2  20 p  +  2*Jqr  sin2  20 p  j 
q  +  r  +  (q-r) (m'm  cos 20 p  +  m'02  sin 20 p ) 


(93) 


and 


"*02  (0P)  = 


(q  -  r)  sin  20 p  +  m'02  ^( q  +  r) sin2  20 p  +  2^Jqr  cos2  20 p  j 
q  +  r  +  (q  —  r)  (m'm  cos  26 p  +  m'02  sin  26 p  ) 


Now,  Equations  (93)  and  (94)  are  plugged  into  Equation  (45).  This  gives 
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(94) 


(< q-r)cos20p  +  m'0l((q  +  r)  cos2  20 p  +  2<s[qr  sin2  20  \ 

1+  Sj  - T - T - 

q  +  r  +  (q  -  r)  ( m'm  cos  20 p  +  m'02  sin  20 p  1 

M»,)  =  L,  (95) 

(g  -  r)  sin  20 p  +  m'02  ( (q  +  r)  sin2  20 p  +  2<Jqr  cos2  20 p  ] 

+s  - - - - 

v  2  g  +  r  +  {q-r ){m'm  cos  20 p  +  m02  sin 20 p )  ^ 

where  the  only  assumption  made  is  that  circular  polarization  in  the  scene  is  negligible  and 
therefore  s3  =  0 .  This  is  a  good  assumption  for  passive  remote  sensing  applications  since 
the  degree  of  circular  polarization  is  typically  very  small  [9].  .v,  and  s2  then  are  solved 
for  using  Equations  (55)  and  (56)  and  normalizing  by  .s0 .  This  gives 


si  = 
1  m 


~(q  -  f)((m'02  - 1  )(q  +  r)sl  +  2m’0lm'02  y[qrs2  ) 

(q  +  r)2  -  Wq2 (q-r)2  +  Am'mqrsx  +  2m'02qsfqrs2  +  2m'02r*jqrs2 


~{q  ~  r )  ( (m02  - 1  ){q  +  r)s2  +  2 m01m02  Jqrsl ) 

s2  = - - - - - .  (97) 

m  (q  +  r)  -m'0l(q-r)  +4m'02qrs2  +  2 m'01qyjqrs1  +2 m'01r^jqrs1 

While  the  form  of  sx  and  s2m  are  complicated,  it  is  interesting  to  see  what  happens  when 

the  Mueller  matrix  of  the  system  (excluding  the  polarizer)  is  assumed  to  be  close  to  or 
exactly  the  identity  matrix.  In  the  case  of  an  identity  system  Mueller  matrix,  m'0l  and 

m'02  are  equal  to  zero.  Using  this  in  Equation  (93),  mm  (0  )  reduces  to 


”hn  (0„) 


(q  -  r)  cos  20 p 


-  D  cos  2 <9 


™m(0n) 


(q  -  r )  sin  20 p 


-  Dsin  20., 
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Using  this  to  solve  for  sx  and  s2m  gives 


5  4,(0)-A„(90)  (1+^)4 -(1-^)4  2s.PL,  sD 

lm  Lm(0)  +  Lm(90)  (\  +  s,D)L+a-s,D)L  2L 

and 

5  4,(45)-Lm(135)  q  +  s2D)L, -(f-sJVL,  sD  (10i) 

2m  Lm(45)  +  Lm(135)  (l  +  s2D)Ls+a-s2D)Ls  2  Ls 

This  means  that  if  the  polarizer  is  non-ideal  and  the  rest  of  the  system  has  no  effect  on  the 
polarimetric  state  of  the  light,  then  the  measured  Stokes  vectors  will  simply  be  the  true 
Stokes  vectors  multiplied  by  the  diattenuation  of  the  polarizer.  This  is  a  good 
approximation  since  the  components  of  the  system  consist  mainly  of  lenses  and  mirrors 
which  have  a  small  polarization  effect.  The  one  element  still  in  question  is  the  beam 
splitter  which  might  have  a  more  significant  effect  since  it  is  more  reflective  and  angled 
with  respect  to  incident  light. 

Simulation  of  Polarimetric  Performance  with  non-ideal  Polarizer 
It  is  hard  to  derive  useful  information  from  Equations  (96)  and  (97),  although  it  is 
interesting  to  note  that  if  the  polarizer  is  ideal,  that  is  q  =  1  and  r  =  0 ,  then  =  .v,  and 

s2m  =  s2  as  expected.  The  easiest  way  to  make  sense  of  Equations  (96)  and  (97)  is  to 
plug  in  realistic  values  for  q  and  r  since  the  extinction  ratio  of  the  polarizer  is  known, 
and  to  plug  in  a  range  of  physically  acceptable  values  for  m'01  and  m'02 .  Then  the  result 

from  various  input  polarization  states  can  be  simulated.  These  results  can  then  be  plotted 
and  compared  to  actual  data  collected  from  the  Hypercam.  Since  the  band-averaged 
extinction  ratio  of  the  instrument  polarizer  is  known  to  be  approximately  400,  let  q  =  0.8 
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and  r  =  0.002 .  This  gives  an  extinction  ratio,  q  /  r  equal  to  400.  Possible  values  for 
m'm  and  m'02  range  from  -1  to  1.  For  the  first  case,  horizontally  polarized  light  will  be 
simulated,  and  therefore,  sx  - 1  and  s2  -  0  is  input  to  the  system.  Doing  this  reduces 
Equations  (96)  and  (97)  to 


i  nn?  0  010 

s.  =1.005 - 

lm  1.010  -m' 


01 


(102) 


and 


s2 


-0.995m'mm'02 


10.025 +  m;i-9.925m;22  ' 
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Figure  49  shows  a  plot  of  sx  in  which  m'0l  was  allowed  to  range  from  -1  to  1.  The  plot 
confirms  that  for  values  of  m'01  close  to  zero,  the  measured  sx  is  approximately  the 
diattenuation  of  the  instrument  polarizer.  A  plot  of  the  measured  s2m  is  more  difficult 
since  it  depends  both  on  m'01  and  m'02 .  To  gain  a  reasonable  picture  of  how  s2  varies 


with  the  values  of  m'm  and  m'02 ,  a  3-D  plot  was  generated  with  both  m'0l  and  m'0 


*02 


ranging  from  -0.7  to  0.7.  As  expected,  the  measured  value  of  s2  is  approximately  zero 


when  m'm  and  m'02  are  close  to  zero. 


The  derivations  thus  far  in  this  section  have  used  the  Pickering  Method  for 
calculation  of  s,  and  sq  .  This  was  done  to  illustrate  the  math  without  having  the 

1  m  zm 

equations  become  overly  long  and  complicated.  By  doing  so,  slm  and  s2m  can  never 

have  a  value  greater  than  one.  It  is  not  hard,  however,  to  extend  the  simulations  to 
incorporate  the  modified  Pickering  Method.  It  will  be  necessary  to  do  so  since  all  of  the 
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results  reported  in  this  thesis  make  use  of  this  method  and  using  this  method  allows  the 
possibility  of  an  s,  or  s2  measurements  greater  than  one  due  to  the  non-ideal  effects 

of  the  polarizer.  Figure  50  shows  the  results  of  the  same  simulation  where  an  Sj  =  1  and 
s2  =0  is  input  into  the  system  and  the  measured  slm  and  s2m  are  determined  using  the 
Modified  Pickering  method.  In  this  case,  the  measured  slm  also  becomes  a  function  of 
both  m'j  and  m'02  and  so  a  3-D  plot  is  used  and  m'm  and  m'2  are  allowed  to  range  from  - 
0.7  to  0.7.  Figure  50  (a)  shows  that  the  measured  sx  has  a  weak  dependence  on  m'02 
and  is  roughly  linearly  dependent  on  m'm  for  reasonable  values  of  m'01  and  m'02 .  Figure 
50  (b)  shows  that  m'm  and  m'02  both  have  a  similar  effect  on  the  measured  s2  ,  but  for 
reasonable  values  of  m'm  and  m'02 ,  s2m  is  approximately  zero  as  expected. 
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Figure  49.  a.)  Simulation  of  the  measured  Sim  Stokes  value  when  horizontally  polarized  light  is 
incident  on  the  system  and  m’0i  is  allowed  to  range  from  -1  to  1.  b.)  Simulation  of  the  measured  s2m 
Stokes  values  when  m’0i  and  m’02  are  allowed  to  range  from  -0.7  to  0.7.  The  extinction  ratio  of  the 

instrument  polarizer  is  set  to  400. 
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Figure  50.  Same  simulation  for  a.)  Sim  and  b.)  S2m  as  in  Figure  49,  except  using  the  Modified 
Pickering  method  for  calculation  of  the  measured  slm  and  s2m.  Notice  that  the  value  of  the  measured 

Si  can  possibly  exceed  one  using  this  method. 

Another  interesting  simulation  to  run  is  for  partially  polarized  light.  If  the 
instrument  polarizer  were  ideal,  the  amount  of  s2  input  would  not  affect  the  measured 
value  for  v,  and  vice-versa.  Since  the  polarizer  is  not  ideal,  however,  the  two 
equations  become  coupled  and  the  s2  input  will  affect  the  measured  value  for  sx  and 
vice-versa.  The  degree  to  which  an  s2  input  affects  the  measured  value  of  s,  is 
dependent  on  the  magnitude  of  m'0l  and  m'02 .  Figure  51  shows  the  simulation  for  the 
measured  value  of  sx  when  the  input  Stokes  values  of  sx  and  s2  are  allowed  to  vary 
from  0  to  1.  Figure  51  (a)  is  for  m'm  =  0.05  and  m'02  =  0.05  and  Figure  51  (b)  is  for 
irim  -  0.5  and  m'02  =  0.5 .  The  results  show  that  the  input  s2  has  very  little  effect  on  the 
measured  sx  when  m'm  =  0.05 ,  but  when  mm  =  0.5  the  input  s2  does  have  a  non¬ 
trivial  effect  on  the  measured  s,  value. 

'm 
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a.) 


Figure  51.  Measured  value  of  slm  when  the  input  Stokes  values  are  allowed  to  range  from  0  to  1  for 
a.)  moi  and  mo2  equal  to  0.05  and  b.)  m0i  and  mo2  equal  to  0.5. 

Finally,  one  last  simulation  is  provided  to  explore  the  dependence  of  the  measured 

Stokes  values  on  the  extinction  ratio  of  the  instrument  polarizer.  To  do  this,  q  and  r  are 

allowed  to  vary  from  0  to  1.  The  input  Stokes  parameters  are  =  1  and  s2  —  0  while 

m'0 j  and  m'02  are  both  set  to  0.05.  Figure  52  shows  the  results.  A  value  of  q  close  to  1 
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and  a  value  of  r  close  to  zero  corresponds  to  a  polarizer  that  is  nearly  ideal.  As  q  moves 


away  from  1  and  r  moves  away  from  0,  the  measured  value  of  sx  decreases  from  the 
true  input  value. 
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Figure  52.  Shows  the  dependence  of  the  measured  Sim  values  on  the  extinction  ratio  of  the  instrument 
polarizer.  The  values  of  q  and  r  are  allowed  to  vary  from  0  to  1.  The  input  Stokes  parameters  are 
Si  =  1  and  s2  =  0  and  m0i  and  moi  are  both  equal  to  0.05. 

The  simulations  shown  above  provide  insight  into  the  effects  of  a  non-ideal 
polarizer  and  a  system  Mueller  matrix  that  is  not  assumed  to  be  the  identity  matrix,  on  the 
quality  of  the  measured  Stokes  parameters.  According  to  Chipman,  in  the  case  where  a 
system  consists  of  weakly  polarizing  elements  such  as  lens  surfaces  and  mirror  surfaces, 
the  polarization  properties  are  not  zero  due  to  the  Fresnel  equations,  antireflection 
coatings,  or  mirrored  surfaces,  but  the  effects  are  often  well  below  5  percent  [14].  For  an 
m'0l  and  m'02  within  ±  5%  of  zero  and  a  horizontal  input  Stokes  vector,  the  values  for 
s.  range  from  0.992  to  0.998  and  v-  ranges  from  -0.0002  to  0.0002.  These  values 
correspond  to  a  ±0.3%  error  in  slm  and  a  ±0.02%  error  in  s2m  .  In  order  to  gain  further 
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insight  into  what  the  m'm  and  m'02  elements  of  the  system  Mueller  matrix  physically 

mean,  the  system  itself  is  treated  as  a  weak  diattenuator  which  means  that  its  Mueller 
matrix  is  given  by  Equation  (92).  Since  the  off-axis  elements  in  Equation  (92)  are 
assumed  to  be  small,  the  first-order  Taylor  series  expansion  can  be  taken.  After 
normalizing  the  Mueller  matrix,  this  gives 

f  i  K  d'45 

M'«  D'»  1  ° 

£>45  0  1 

[d'r  0  0 

where  m'01  ~  D'H  determines  the  transmission  of  a  horizontal/vertical  input  state  and 
mo2  ~  £>45  determines  the  transmission  of  a  457135°  input  state.  In  other  words,  the  m'01 
and  m'02  elements  determine  the  diattenuation  of  the  HSI  instrument  itself. 
Determination  of  the  HSI  System  Diattenuation 

Now,  using  the  actual  measured  response  of  the  system  for  the  12  different  input 
Stokes  vectors,  the  system  model  can  be  used  to  solve  for  the  m'01  and  m'02  elements  of 

the  HSI  system's  Mueller  matrix.  This  will  then  allow  for  the  calculation  of  the  HSI 
system's  diattenuation.  As  a  reminder,  when  referring  to  the  HSI  system,  the  polarizer  is 
being  excluded.  Equations  (96)  and  (97)  are  re-stated  here  for  reference: 

-(q  -  r ) ( (m'0 2  -  \)(q  +  r)sx  +  2 m'mm'02  y [qrs2 ) 

^1  =  - - - - - - - - - - - 

lm  ( q  +  r )2  - m'^fq - r)2  + 4 m'01qrs1  +2m'02q<Jqrs2  +2 m'02r^Jqrs2 

~{q  ~  r )  ( (m'0 2  - 1  )(q  +  r)s2+  2 m'mm'02  <J, qrs x ) 

^2  =  - - - .  - - - ^  - 

m  (q  +  r)2  -m'^2(q-r)2  +4m'02qrs2  +2m'0lqsfqrsl  +2m'mr^jqrsl 


o' 

0 

1 


'  A 


(104) 
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where  sx  and  s2m  are  the  measured  Stokes  parameters  for  a  given  input  sx  and  s2 .  Also 

recall  that  these  equations  were  derived  using  the  Pickering  method  whereas  the  same 
equations  derived  using  the  Modified  Pickering  method  are  similar  but  longer  and  more 
complicated.  The  results  described  in  this  section  use  the  Modified  Pickering  method 
equations.  To  find  the  best  solution  for  m'01  and  m'02 ,  the  Nelder-Mead  simplex  direct 

search  algorithm  was  implemented  to  minimize  the  error  between  the  true  measured 
values  for  sx  and  s2  and  the  values  generated  with  above  equations  while  allowing 

m'0l  and  m'02  to  range  from  -1  to  1.  The  results  give  a  value  of  m'm  =  0.14  and 

mo2  =  -0.41 .  According  to  Chipman  [14],  a  diattenuation  parameter  set,  d ,  can  be 

defined  for  linearly  polarized  light  as 

jmn,  mm )  ,  ,  ,  , 

d  = - =  {m'm  m'02}  =  [dH  d45],  (105) 

moo 

which  characterizes  the  dependence  of  the  transmission  on  the  incident  polarization  state. 
The  linear  diattenuation  of  any  Mueller  matrix  is  then  given  by 

D(M)  =  =  /m't  +  m'2  ,  (106) 

T  +  T  1 

max  min 

Utilizing  this  to  determine  the  linear  diattenuation  of  the  HSI  system  based  on  the 
calculated  values  of  m'0l  and  m'02  gives 


D(MHS! )  =  ^(0.14)2+(-0.41)2  =  0.43 . 


(107) 


The  extinction  ratio  of  the  HSI  system  can  also  be  determined  with 

1 +  D  1  +  .0.43 


ER  = 


1  -D  1-0.43 


=  2.51. 


(108) 
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To  visualize  how  well  the  calculated  values  of  s,  and  .s'  match  the  measured  values, 

'm  zra 


Figure  53  plots  them  together. 
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Figure  53.  The  measured  values  of  Sim  and  s2m  for  the  12  input  Stokes  vectors  are  plotted  with  circles 
and  the  calculated  values  for  Sim  and  s2mare  plotted  with  squares.  The  input  curve  is  plotted  as  a 

dashed  line. 
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VI.  Conclusions 


Polarimetric-hyperspectral  remote  sensing  is  a  promising  field  that  brings  two 
traditionally  independent  modalities  together  to  enhance  scene  characterization 
capabilities.  The  Telops  P-HSI  instrument  combines  these  technologies  in  a  way  that  has 
never  been  done  before  and  provides  a  combined  imaging  and  spectral  capability  that  is 
considered  state-of-the  art.  The  Defense  Threat  Reduction  Agency  (DTRA)  funded 
research  at  AFIT  to  leverage  this  capability  to  better  inform  radiation  transport  models 
that  will  allow  for  more  effective  location  of  ionizing  radiation  sources  within  the  scene. 

The  desire  for  highly  accurate  data  requires  careful  calibration  of  the  instrument. 
In  this  thesis,  a  mathematical  calibration  framework  was  developed  that  links  standard 
Fourier  transform  spectrometer  (FTS)  calibration  with  standard  polarimetric  calibration 
in  a  simple,  straightforward  manner.  Building  this  framework  was  key  for  two  reasons. 
First,  it  was  essential  for  ensuring  the  calibration  methodology  was  sound  and  it  provided 
a  framework  for  understanding  the  influence  of  various  instrument  parameters  (both  ideal 
and  non-ideal)  on  ultimate  calibration  performance.  Second,  this  development  makes  it 
easy  for  both  FTS  and  polarimetric  experts  to  understand,  and  paves  the  way  for  use  of 
this  method  for  similar  instruments  in  the  future. 

Due  to  the  complexities  of  the  system  and  operation  in  LWIR  waveband,  three 
different  types  of  calibration  are  required.  These  are:  radiometric  calibration,  off-axis 
spectral  calibration,  and  polarimetric  calibration.  As  a  result  of  this  thesis,  all  three 
calibrations  are  now  successfully  accomplished.  This  means  that  highly  accurate 
radiometric,  spectral,  and  polarimetric  data  can  be  generated  with  the  Telops  instrument. 
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This  will  enable  more  accurate  scene  characterization,  which  will  provide  more  accurate 
information  to  the  radiation  models  developed  under  DTRA  sponsorship. 

While  all  three  types  of  calibration  were  demonstrated  in  this  thesis,  the  primary 
focus  was  on  polarimetric  calibration.  This  was  the  first  time  the  Telops  instrument  has 
been  used  with  the  additional  polarimetric  capability,  and  therefore,  this  is  the  first  time 
polarimetric  calibration  has  been  performed  on  the  instrument.  The  polarimetric 
calibration  effort  showed  that  by  taking  calibration  blackbody  data  at  each  of  the  four 
polarizer  angles,  and  solving  for  the  spectro-polarimetric  gain  and  offset  at  each  angle, 
the  full  polarimetric  calibration  is  achieved  with  the  only  assumption  being  that  the 
polarizer  is  ideal,  i.e.,  it  has  a  diattenuation  of  one.  The  polarimetric  calibration  was  then 
expanded  to  include  the  more  general  case  of  a  non-ideal  polarizer  and  the  methodology 
necessary  to  correct  for  and  quantify  the  non-ideality  was  developed.  This  involved 
experimental  determination  of  the  system  matrix,  data  reduction  matrix,  and  Mueller 
deviation  matrix  for  each  pixel  at  each  waveband. 

The  polarimetric  calibration  was  then  performed,  first  assuming  an  ideal  polarizer, 

and  then  making  the  additional  correction  for  a  non-ideal  polarizer.  The  additional 

correction  for  a  non-ideal  polarizer  was  difficult  considering  the  fact  that  the 

diattenuation  of  the  instrument  polarizer  is  approximately  0.995,  and  therefore,  the 

necessary  correction  is  only  roughly  0.5%.  The  results  showed  that  the  correction  did 

improve  the  error  by  0.56%  or  less  on  average.  The  results  also  showed,  however,  that 

there  was  some  residual  error  left  over  after  the  correction  and  that  the  correction  is 

ultimately  limited  by  the  noise  in  the  measurements,  the  non-linearity  of  the  detector,  and 

the  accuracy  of  the  experimental  setup.  Each  of  these  limitations  were  discussed  and  the 
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noise  statistics  of  the  measurements  were  explored.  Since  the  correction  is  required  on  a 
pixel-by-pixel  and  band-by-band  basis,  the  noise  was  determined  on  a  pixel-by-pixel  and 
band-by-band  basis.  It  was  shown  that  some  of  the  residual  error  in  the  polarized  data 
can  be  attributed  to  noise  in  the  measurements.  It  was  also  shown  that  a  primary  source 
of  the  residual  error  in  the  data  was  likely  due  to  the  thermal  stability  of  the  6-inch 
blackbody  source.  Additionally,  it  was  shown  that  there  is  currently  a  fundamental 
limitation  on  all  data  sets  which  is  caused  by  the  non-linearity  of  the  detector.  This 
causes  a  checker  board  pattern  on  all  images  with  low  signal  levels.  This  same  checker 
board  pattern  shows  up  in  the  Mueller  deviation  matrix  which  represents  the  non¬ 
idealities  of  the  system.  Determining  the  exact  cause  and  proper  correction  for  this 
checker  board  pattern  is  a  subject  of  on-going  research. 

In  the  process  of  deriving  the  full  spectro-polarimetric  calibration  methodology 
for  the  instrument,  a  polarimetric  system  model  was  developed.  This  system  model  maps 
an  input  scene  Stokes  vector  to  the  Stokes  vector  that  is  actually  measured  by  the 
instrument  as  a  function  of  the  extinction  ratio  of  the  instrument  polarizer,  the 
diattenuation  of  the  HSI  system  (not  including  the  polarizer),  and  the  input  stokes 
parameters.  This  allowed  the  system  performance  to  be  simulated  as  a  function  of  the 
polarimetric  system  parameters  and  gave  insight  into  how  each  parameter  affects  the 
accuracy  of  measured  data.  It  was  shown  that  in  the  case  of  a  small  HSI  system 
diattenuation,  the  measured  data  was  accurate  and  the  primary  limitation  was  imposed  by 
the  polarizer  extinction  ratio.  A  large  system  diattenuation,  however,  imposed  error  on 
the  measured  data  which  must  be  corrected  with  the  calibration  methods  described  in  this 


thesis.  The  actual  diattenuation  of  the  HSI  system  (not  including  the  polarizer)  was  then 
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determined  using  the  system  model  and  the  measured  data  from  the  experiment.  The 
diattenuation  of  the  system  was  found  to  be  0.43  which  corresponds  to  an  extinction  ratio 
of  2.51.  This  is  small  compared  with  the  polarizer's  extinction  ratio  of  400;  however,  it 
must  be  calibrated  out  to  achieve  the  most  accurate  data  from  the  instrument.  This  thesis 
demonstrated  the  methodology  required  to  calibrate  this  effect  out. 

In  addition  to  calibration,  several  figures  of  merit  were  determined  for  the  system. 

These  included  the  NESR  and  NESDoLP  for  both  unpolarized  and  polarized  light.  For 

unpolarized  light,  the  NESR  was  reported  at  four  instrument  polarizer  angles  (0°,  45°, 

90°,  135°).  The  results  showed  that  the  NESR  at  45°  and  135°  were  equal  across  the 

band,  whereas  the  NESR  for  0°  was  lower  and  the  NESR  for  90°  was  higher.  The  values 

reported  are  for  a  spectral  resolution  of  16  cm'1,  viewing  a  55  °C  blackbody  source  with  a 

200  jus  integration  time.  The  band-averaged  values  were  30.6  nW/(cm2sr)  cm"1,  33.6 

nW/(cm2sr)  cm"1,  33.6  nW/(cm2sr)  cm'1,  and  36.9  nW/(cm2sr)  cm'1  for  the  0°,  45°,  135°, 

and  90°  measurements  respectively.  This  result  verified  that  the  system  was  operating 

within  the  specifications  provided  by  the  manufacturer.  The  fact  that  the  NESR  for  the  0° 

and  90°  measurements  were  split  while  the  NESR  for  the  45°  and  135°  were  equal  is 

indicative  of  the  diattenuation  of  the  HSI  system  after  the  light  passes  through  the 

polarizer.  The  band-averaged  NESDoLP  for  the  unpolarized  source  was  0.39%.  The 

actual  calculated  DoLP  for  the  unpolarized  source  was  0.51%.  This  indicates  that  a 

residual  error  or  bias  of  0.12%  remained  in  the  data  after  calibration.  This  residual  error 

was  attributed  to  either  thermal  drift  in  the  reference  source  or  detector  non-linearity. 

The  band-averaged  NESR  for  the  polarized  source  was  29.8  nW/(cm2sr)  cm'1,  31.8 

nW/(cm2sr),  39.1  nW/(cm2sr),  and  41.0  nW/(cm2sr)  for  the  45°,  90°,  0°,  and  135° 
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measurements  respectively.  This  again  was  within  the  specifications  provided  by  the 
manufacturer.  The  measured  DoLP  for  the  polarized  source  was  100.48%  before 
correction  for  the  non-ideal  polarizer  and  99.92%  after  correction.  The  DoLP  above 
100%  was  attributed  to  noise  in  the  measurements  as  well  as  the  possible  thermal  drift  of 
the  reference  blackbody  source.  In  this  case,  the  correction  improved  the  accuracy  of  the 
DoLP  measurement  by  0.56%,  but  left  a  residual  error  of  0.62%. 

Ultimately,  this  thesis  provides  solid  evidence  that  the  Telops  instrument  is 
performing  as  expected  from  a  polarimetric  standpoint  and  this  is  something  that  can  be 
presented  to  the  polarimetric  community  for  critique.  This  will  then  allow  future 
researchers  to  have  confidence  in  their  results  while  focusing  more  on  the  applications  of 
the  data  such  as  material  identification,  target  detection,  anomaly  detection,  disturbed 
earth  detection,  and  3-D  scene  reconstruction.  All  of  which  are  important  capabilities 
that  the  DoD  and  IC  wish  to  develop  and  improve  upon.  With  a  few  simple 
improvements  to  the  experimental  setup,  the  residual  errors  reported  in  this  thesis  should 
be  reduced  to  negligible  values.  The  next  section  discusses  recommendations  to  improve 
the  setup  for  future  calibration  efforts  in  order  to  achieve  the  best  results. 

Future  Recommendations 

Three  iterations  of  the  calibration  experiment  were  completed,  each  one 
improving  results  over  its  predecessor.  The  final  experiment  was  done  in  a  way  that 
minimized  errors  due  to  the  setup  as  much  as  possible  within  the  limitations  of  the 
available  equipment  at  the  time.  After  a  thorough  analysis  of  the  data  and  a  better 
understanding  of  the  sources  of  residual  error,  there  are  still  several  measures  that  can  be 
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taken  to  improve  the  accuracy  of  the  calibration.  Since  the  error  introduced  by  assuming 
an  ideal  instrument  polarizer  is  so  small,  the  experimental  setup  must  achieve  extremely 
high  accuracy  in  order  to  provide  the  additional  correction  for  the  non-ideal  nature  of  the 
polarizer.  Here  are  several  recommendations  for  achieving  higher  accuracy  with  the 
experimental  setup  and  thus  better  results  from  the  calibration. 

First,  it  is  necessary  to  know  the  input  Stokes  vectors  with  very  high  accuracy.  If 
a  linear  wire-grid  polarizer  is  used  as  in  the  above  experiment,  this  means  careful 
characterization  of  the  spectral  diattenuation  of  the  polarizer  and  high  confidence  that 
there  are  no  blemishes  in  the  surface  of  the  polarizer  that  affect  it  spatially.  For  this 
thesis,  the  manufacturer's  data  sheet  was  used  for  determination  of  the  spectral 
diattenuation;  however,  this  characterization  could  be  done  experimentally,  right  before 
the  calibration.  Additionally,  use  of  an  accurate,  motor  driven  rotation  stage  would  be 
useful  for  decreasing  the  uncertainty  in  the  external  polarizer  orientation  and  increasing 
the  accuracy  of  the  input  Stokes  vectors. 

Second,  the  6-inch  blackbody  should  be  eliminated  from  the  setup  due  to  its  poor 
thermal  stability  performance  ( ±  0.5  °C).  This  is  an  unacceptable  thermal  drift,  which  is 
likely  the  primary  source  of  the  residual  error  seen  in  the  data.  The  6-inch  blackbody  was 
used  to  ensure  that  the  entire  FOV  of  the  polarizer  was  filled  by  the  blackbody,  however, 
with  careful  setup  procedures,  the  6-inch  blackbody  could  be  replaced  with  a  more  stable 
2-inch  blackbody  ( ±  0.003  °C). 

Third,  to  improve  the  noise  statistics,  it  would  have  been  good  to  increase  the 
difference  between  the  hot  and  cold  data  sets.  That  way,  when  they  are  differenced,  the 


SNR  is  better.  This  combined  with  a  larger  number  of  frame  averages  would  improve  the 
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calibration.  Since  the  system  matrix  was  shown  to  be  flat  spectrally,  a  lower  spectral 
resolution  could  be  used  to  help  increase  SNR  as  well. 

Fourth,  the  external  polarizer  could  be  moved  much  closer  to  the  Telops 
instrument.  In  this  thesis,  the  external  polarizer  was  placed  at  a  range  that  would  allow 
the  polarizer  to  be  in  focus.  This  isn't  necessary  however,  and  by  moving  it  closer  to  the 
instrument,  there  is  a  smaller  possibility  of  the  target  blackbodies  causing  thermal  effects 
on  the  external  polarizer.  Moving  it  closer  would  also  allow  more  pixels  on  the  FPA  to 
be  calibrated  at  one  time. 

Finally,  a  method  should  be  developed  to  perform  the  polarimetric  calibration  on 
the  entire  FPA.  This  capability  is  limited  by  the  small  extent  of  the  external  polarizer, 
which  prevents  it  from  filling  the  entire  FOV  of  the  system.  One  suggestion  would  be  to 
perform  the  calibration  on  one  region  of  the  FPA  at  a  time  until  the  entire  FPA  is 
mapped.  The  system  matrix  is  not  expected  vary  much  in  the  spatial  dimension  and  so  it 
should  not  be  hard  to  map  the  entire  FPA  with  this  method. 
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